Spin-charge separation in one-dimensional fermion systems 
beyond the Luttinger liquid theory 
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We develop a nonperturbative zero-temperature theory for the dynamic response functions of 
interacting one-dimensional spin-1/2 fermions. In contrast to the conventional Luttinger liquid 
theory, we take into account the nonlinearity of the fermion dispersion exactly. We calculate the 
power-law singularities of the spectral function and the charge and spin density structure factors 
for arbitrary momenta and interaction strengths. The exponents characterizing the singularities 
are functions of momenta and differ significantly from the predictions of the linear Luttinger liquid 
theory. We generalize the notion of the spin-charge separation to the nonlinear spectrum. This 
generalization leads to phenomenological relations between threshold exponents and the threshold 
energy. 
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I. INTRODUCTION 

One of the aspirations of condensed matter physics 
is to understand the physical properties of interacting 
many-body systems. A successful example is Landau's 
Fermi liquid theory^ which provides a comprehensive 
framework for the description of repulsively interact- 
ing fermions. This theory is based on the concept of 
fermionic quasiparticles, which can be regarded as phys- 
ical fermions surrounded by a cloud of particle-hole ex- 
citations. The quasiparticles carry the same quantum 
numbers as the physical fermions but generally have a 
different effective mass. Even if the interactions between 
the physical fermions are strong, phase space constraints 
near the Fermi surface limit the scattering rates of quasi- 
particles and they become stable towards low energies. 
In this limit, the system can be described by a theory of 
nonintcracting fermionic quasiparticles. 

Let us compare some properties of nonintcracting 
Fermi gases and interacting Fermi liquids in three di- 
mensions. In a Fermi gas, the occupation number n(k) = 
(iplipk) which measures the number of physical fermions 
with momentum k, jumps from one to zero at the Fermi 
surface \k\ — kp. For a Fermi liquid the amplitude of the 
discontinuity at \k\ = kp is reduced to a positive value 
Z < 1, but remains nonzero. The so-called quasiparticle 
residue Z is a measure of the overlap between the physi- 
cal fermions and the quasiparticles. Another example is 
the spectral function A(k,ui), defined below in Eq. ©. 
For a Fermi gas with spectrum e(/c), the spectral function 
is A(k, ui) = 6 [ui — e(k)] . For a Fermi liquid, the <5-peak at 
energies close to the Fermi level evolves into a symmetric 
Lorentzian peak with a width proportional to (\k \ — kp) 2 . 

The physical properties of interacting systems are dras- 
tically different in one dimension. The effect of inter- 
actions is nonperturbative. For example, an arbitrarily 
weak repulsion in one dimension leads to zero residue, 
Z = 0, suggesting an absence of fermionic excitations. 



Instead, the elementary excitations are thought to be 
better represented by quantized waves of density obeying 
Bose statistics^ For low energies, the spectrum e{k) of 
the physical fermions can be linearized around the two 
Fermi points ±kp. The interacting fermionic theory can 
then be mapped onto a theory of noninteracting bosons^ 
and all correlation functions can be calculated exactly^ 
The universality class formed by gaplcss interacting one- 
dimensional systems at low energies is called Luttinger 
liquid (LL)£ It is entirely characterized by the velocity 
of density waves v and the Luttinger parameter K which 
depends on the interaction strength^ The distinct nature 
of these systems compared to their higher-dimensional 
counterparts manifests itself in many observables. The 
spectral function A(k,Lu), for example, instead of becom- 
ing a Lorentzian, displays asymmetric power-law diver- 
gencies. 

The contrast between one-dimensional and higher- 
dimensional interacting systems becomes even stronger 
for spinful systems. The eigenmodes of a spinful LL 
are spin-carrying and charge-carrying density waves with 
linear dispersion^ In general, the velocity of these two 
types of excitations, v s and v c respectively, can be very 
different. The introduction of a physical fermion with 
charge and spin into the liquid leads to the formation 
of independent spin and charge density waves. This 
phenomenon is called spin-charge separation and can be 
probed in experiments^— When turning on the interac- 
tions, the 5-function singularity in the spectral function 
of a noninteracting system, A(q,u) = S(lo — vfq), splits 
into two power-law singularities at the threshold ener- 
gies for spin and charge density waves ; 12 ' 13 A(q, to) cx 
(a; — v s q)~^ (uj — v c q)~' tJ ' with exponents fj. c,s depending 
only on the Luttinger parameter. 

The linearization of the generic spectrum of particles 
is the crucial simplification leading to the LL theory and 
it has a substantial impact on density excitations. For an 
LL, the charge and spin density structure factors defined 
below in Eq. (JTJ) have sharp peaks, S(q, co) cx 5(lj — v c \q\) 
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and S zz (q,uj) oc S(u — v s \q\) respectively, because a linear 
spectrum entails a one-to-one correspondence between 
the energy and the momentum of density excitations. 
This is the reason why the eigenmodes of an LL are den- 
sity excitations. Also note that spin and charge modes 
are entirely decoupled. 

For a quadratic spectrum, in contrast, a density exci- 
tation of fixed momentum q may have energies in a range 
of width proportional to q 2 . The peak in the structure 
factor remains narrow in the sense that q 2 /(v c ^ s q) cx q, 
as long as q remains small compared to kp. This agrees 
with the general idea of irrelevance of the curvature as 
a perturbation to the LL theory^ However, the form 
of the peak is far from being a simple broadening to a 
Lorentzianji 4 - Moreover, a spectrum curvature leads to a 
coupling between the spin and charge modes.— ~— The 
charge density structure factor, for example, acquires a 
peak at energies uj ~ v s q characteristic for the spinon 
excitations! 18 ' 19 

Along with the structure factors, the spectral func- 
tion is also affected by the nonlincarity of the dispersion 
rclationj22r— The energy domains in which the dynamic 
response functions deviate from the predictions of the LL 
theory expand with the particles' momenta tuned away 
from the Fermi points. It was shown, however, that even 
far away from the Fermi points there is a certain uni- 
versality in the behavior of the response functions . 21 ! 22 
The dynamic response functions at arbitrary wave vec- 
tors of single-species fermions or bosons were studied in 
detaiL 14 i 23 " — The case of spin-1/2 fermions, despite its 
practical importance, is studied far less.— ^ 

The experiments on momentum-resolved tunneling be- 
tween ID systems^ 9 - or between ID and 2D systems^ 
as well as angle-resolved photocmission spectroscopy on 
quasi-lD systems , 10 ' 30 " — provide a tool to measure the 
electron spectral function. Moreover, the density struc- 
ture factor can be measured using the Coulomb dra g 23 ' 34 
or neutron scattering on spin chains.— Some of these ex- 
periments may be interpreted, in a limited domain of 
wave vectors, in terms of spin-charge separation! 10 i 11 i 31 
Moreover, experiments using ultracold gases have been 
propose d 36 ' 37 which could allow the observation of spin- 
charge separation in real space. The numerical calcu- 
lation of dynamic response functions has become possi- 
ble using time-dependent density matrix renormalization 
group techniques.— ~— None of the developed methods is 
limited to low energies and momenta, thus prompting the 
question of the nonlinear dispersion effects. 

The main goal of this article is to present quantita- 
tive results for the dynamic response functions (spectral 
function and the density structure factors) for spinful ID 
Fermi systems at arbitrary interaction strength. We ex- 
tend the LL theory by taking into account the nonlinear- 
ity of the fermion spectrum exactly, and we obtain re- 
sults for the dynamic response functions which are valid 
for arbitrary momenta. Moreover, we shall elucidate the 
fate of the spin-charge separation away from the Fermi 
points. Our results apply to a wide range of systems with 



gaplcss spectrum and will be used to track the evolution 
of the dynamic response functions all the way from the 
nonintcracting to the strongly interacting limit. 

We restrict our analysis to spin-1/2 systems at zero 
magnetic field, i.e., without spin polarization. Spin rota- 
tion symmetry then entails SU (2)-invariance of the spin 
degrees of freedom. On the other hand, spin-polarized 
one-dimensional systems are interesting in their own right 
and have been investigated theoretically 4 ^ as well as 
experimentally— A nonzero magnetic field breaks the 
5f7(2)-symmetry and the Zeeman shift leads to different 
Fermi wavevectors for spin-up and spin-down fermions. 
This splits the peaks in the dynamic response functions, 
and for large magnetic fields spin and charge degrees 
of freedom become coupled even within the linear LL 
theory4£ These complications do not arise in the absence 
of a magnetic field. 

In order to calculate the dynamic response functions, 
it is convenient to translate the bosonic spin and charge 
modes into fermionic quasiparticles, spinons and holons. 
For a linear spectrum, the bosonic or fermionic languages 
may be used equally comfortably and both offer their par- 
ticular benefits. The advantage of the former is the di- 
rect relation between the bosonic modes and the density 
response functions. On the other hand, the fermionic 
description connects to the well-known physics of the 
Fermi edge problem! 4 ^— For a nonlinear spectrum, the 
fermionic basis is superior because it avoids divergencies 
arising in the bosonic perturbation theory.— It leads to a 
generalization of the quantum impurity model which was 
used previously to calculate dynamic response functions 
for spinless systems! 21 ' 22 

This article is organized as follows: in Sec. UJ we 
present an overview of our results for the dynamic re- 
sponse functions at zero temperature and the spin-charge 
separation, and point out qualitatively the main differ- 
ences between interacting systems with linear and non- 
linear spectrum. In Sec. IIII| we rephrase the LL theory 
in the basis of fermionic spin and charge quasiparticles 
and reproduce the known results for the spectral func- 
tion. We also discuss how the Hamiltonian changes in the 
presence of a nonlinear spectrum. In Sec. IIV1 we present 
in detail the method for the calculation of threshold sin- 
gularities of dynamic response functions for a nonlincar 
spectrum. We express the threshold exponents of the 
dynamic response functions in terms of scattering phase 
shifts and we calculate the spectral function at its edge 
of support for \k\ — > kp. In Sec. |Vj we construct phe- 
nomenological relations between these scattering phase 
shifts and the shape of the spinon spectrum e s (k) which 
are valid for arbitrary momenta in Galilean invariant sys- 
tem. In Sees. IV Bl and IV C| we use these relations to 
calculate the spectral function A{k, w) and density struc- 
ture factors S zz (k,uj), S h (fc,w) and S(k,w) near their 
respective edges of support for arbitrary momenta. In 
Sec. IVII we construct phenomenological relations fixing 
the exponents of correlation functions near the holon 
spectrum e c (k) in terms of its shape. In Sec. I VIII we 
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apply our general theory to the limits of very strong and 
very weak interactions. For strong interactions, we repro- 
duce the known results for the Hubbard model with in- 
finite interaction. For weak interactions, we complement 
the phenomenological result by a perturbative calculation 
in the basis of free fermions in order to obtain A(k, uS) 
away from the spinon and holon spectra. In Sec. IVI1 we 
estimate the width of the peak of A(k, to) at the holon 
mass shell by investigating the decay rate of holons due 
to their interaction with spinons. Finally, in Sec. I FX! we 
present our conclusions. 



II. QUALITATIVE PICTURE AND RESULTS 

The LL theory is widely used to describe the low- 
energy properties of gapless one-dimensional interacting 
fcrmionic systems £ The restriction to low energies usu- 
ally justifies a linearization of the spectrum of the phys- 
ical fermions around the right and left Fermi points, 
e(k) rj vp(±k — kp). Within this approximation, the 
system remains exactly solvable even for nonzero inter- 
actions and can be cast into a linear theory of nonintcr- 
acting bosonic fields.— These eigenmodes are collective 
density waves which correspond to many-particle excita- 
tions when expressed in terms of the physical fermions. 

One of the notable features of spinful interacting sys- 
tems is the spin-charge separation^ The Hamiltonian of 
the interacting system splits into a sum of two commut- 
ing quadratic terms which act on different Hilbcrt spaces 
and describe the charge and spin degrees of freedom sep- 
arately. For nonzero interactions, the velocities of these 
two types of excitations are different. The injection or 
extraction of a physical particle which carries both spin 
and charge thus leads to the formation of spin and charge 
density waves which separate in space. 

Both the collective nature of the eigenmodes and the 
spin-charge separation are clearly observable in various 
dynamic response functions. The charge and spin density 
structure factors are defined as 

S(k,u) = J dxdte l ^ lkx (p c (x,t)p e (0,0)), 

S~ + (k,u) = J dxdte iuJt ~ ikx (S~(x,t)S + (0,0)) , 

S zz (k,uj) = [ dxdte lut ~ lkx (S z (x,t)S z (0,0)) , (1) 



where p c (x) and S(x) denote the charge and spin density, 
respectively, and S = S x ± iS v . The density structure 
factors measure the linear response of the system at mo- 
mentum k and energy uj to a perturbation which couples 
to the charge or spin density. The charge density pertur- 
bation, for instance, can be created by the absorption of 
a photon. For an LL, one finds S(k,ui) oc 5(u) — v c \k\) 
and S zz (k,uj) = ±S~+(k,uj) cx 8(u - v s \k\). The Dirac- 
S shape of these functions reflects the fact that charge 
and spin density waves are eigenmodes and thus have a 



sharp energy for a given momentum. Note that this is 
a consequence of the linearized spectrum of the physi- 
cal fermions. Moreover, the functions demonstrate that 
charge and spin density waves propagate with velocities 
v c and v s , respectively, which depend on the details of 
the interaction between the physical fermions. 

The spectral function is defined in terms of the re- 
tarded Green's function by A(k,uj) = -ilmG ret (fc,w) 
and can be written as 



roc /*oo 

A(k, u)) = — Re / dt dxe 
n Jo J-oo 

x <{W(M),V>1(0,0)}>. 



iujt—ikx 



(2) 



The operator ipl-(x) creates a physical fcrmion of spin 
(T =t,|- In the absence of a magnetic field, SU(2)- 
symmetry ensures that A(k,u>) is independent of a. 
A(k, uj) represents a different type of dynamic response 
function which measures the response of the system to 
the addition of a physical particle or hole with mo- 
mentum k and energy u>. This function can be deter- 
mined experimentally, for instance, by measuring the 
momentum-resolved tunneling into LLs&ii or by photoc- 
mission spectroscopy ! 10 ' 32 For a nonintcracting system, 
A(k,u>) = 5[oj — e(fe)] defines the spectrum for single- 
particle excitations, e(k). In the presence of interactions, 
customarily described by an LL, this function develops 
power-law singularities at the eigenenergies of spin and 
charge density waves, 



A[k,u> 



s(k)] OC [u) 



.(*)]■ 



»(*) 



(3) 



The qualitative shape of A(k,ui) is shown in Fig. [T] For 
an LL, the spectrum of left- and right-moving spin and 
charge density waves is linear, e ClS (fc) = w c s (±fc — kp), 
and the exponents p°' s are ^-independent. They only de- 
pend on the Luttingcr parameter K c which encodes the 
interaction strength: for a nonintcracting system K c = 1 , 
while repulsive interactions lead to < K c < 1. For 
u) = v c (k — kp), the incoming particle leads to the for- 
mation of a charge density wave with energy uj while 
the spin density wave carries no energy. Similarly, for 
lu = v s (k — kp), the final state contains a spin density 
wave with energy uj and a charge density wave of zero 
energy. Away from these thresholds, the final state may 
contain multiple excitations of nonzero energy in the spin 
and charge sectors. 

In the following, we shall refer to this conventional de- 
scription as the linear LL theory in order to emphasize 
the distinction to the case of nonlinear spectrum. Away 
from the Fermi points, the curvature of the physical spec- 
trum e(k) can no longer be neglected. It is convenient 
to refermionize the system and to express the spin and 
charge density waves in terms of fermionic quasiparticles, 
spinons and holons. Within the linear LL theory, the 
spinon and holon spectra e s (k) and e c (fc) are linear and 
the quasiparticles are nonintcracting. In contrast, for 
nonlinear e(k) spinon and holon spectra become them- 
selves nonlinear and interactions among spinons, among 
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FIG. 1: (Color online) Spinon and holon spectra, e 3 (k) and 
e c (fc), and spectral function A(k, oj) for momenta k > kp- For 
repulsive interactions, a nonlinear fermion spectrum reduces 
E s (k) and increases e c (k) compared to the linear Luttinger 
liquid spectra. The spectral function has a power-law singu- 
larity at the spinon mass shell oj ~ e s (fc) with exponent fi" 
and e a (k) is the edge of support of A(k,uj). The singularity 
at the holon mass shell to « e c (k) is smeared out (dotted line) 
away from kF for nonintegrable systems due to holon decay. 



holons and between spinons and holons come into exis- 
tence. 

Throughout this article, we shall focus on the case of 
zero temperature and repulsive interactions between the 
physical fermions. Repulsion leads to v s < v c , as known 
from the linear LL theory. Let us discuss A{k > kF,ui) 
for nonlinear spectrum. For momenta close to the Fermi 
points, e s (k) bends downwards away from the linear spec- 
trum as depicted in Fig. [T] The spinon spectrum e s (k) 
becomes the edge of support of the spectral function, i.e., 
A(k,uj) = for — e s (k) < oj < e s (k). Near the edge, for 
ui « e s (fc), A(k,uj) is given by Eq. ([3]) with a fc-dependent 
exponent fi s (k). This exponent coincides with the LL 
prediction in the limit k — > kF- In order to calculate the 
spectral function away from kp, we derive phenomcnolog- 
ical expressions which yield the threshold exponent /I s (k) 
at arbitrary k in terms of the properties of e s (k). As the 
latter can be measured or calculated for many systems, 
this provides a useful relation between two independently 
measurable quantities. 

On the other hand, as shown in Fig. [TJ the holon spec- 
trum e c (fc) bends upwards away from kF- For generic 
nonintegrable systems, the singularity of the spectral 
function at oj w e c (fc) becomes smeared out because 
quasiparticle interactions allow a decay of holons. There- 
fore, the power-law behavior at the holon mass shell of 
the form © manifests itself only at k — > kp. It turns 
out that for nonlinear e c (fc), the exponent /i c is different 
form the LL prediction even at k — > kp. The reason is 
that the leading quadratic curvature of e c (k) introduces 
a new energy scale (k — kp) 2 /(2m*), where m* is the 
effective mass. The modified exponent only holds in an 
energy window of this width around e c (k). Since for a 
strictly linear spectrum, the width shrinks to zero, this 
does not contradict the linear LL theory. Thus for the 
charge mode the effect of the spectrum nonlinearity is 
similar to the spinless fermion case^ 

The strict spin-charge separation of the linear LL the- 
ory no longer exists once the band curvature is taken 



into account. The Hamiltonian does not consist of com- 
muting spin and charge terms any more. The dynamic 
structure factors for the charge and spin density, S(k,oj) 
and S zz (k,ui) respectively, cease to be 5-functions. Nev- 
ertheless, spin-charge separation continues to hold in a 
weaker sense: for arbitrary momentum k, the power-law 
singularity of the spectral function at its edge of support 
is determined by states where a spinon carries the entire 
energy e s (k) and has the velocity de s (k)/dk. The holon, 
in contrast, has zero energy and a strictly higher veloc- 
ity v c . Slightly away from the edge, the excess energy is 
used to generate additional low-energy particle-hole pairs 
in the holon sector, but no additional spinons. Moreover, 
spinon decay (in contrast to holon decay) remains forbid- 
den by energy and momentum conservation. Therefore, 
the injection or extraction of a particle with momentum 
k and energy near the threshold, oj w e s (fc), still forms 
spatially separating spin and charge density waves as in 
the linear LL theory. However, in contrast to the linear 
LL theory, this is no longer true for energies far away 
from the threshold. 

The results for the spectral function and the density 
structure factors near the edges of support for arbitrary 
momenta are summarized in Table HI 

III. REFERMIONIZATION OF THE 
LUTTINGER LIQUID THEORY 

Let us start by recapitulating the results for the spec- 
tral function of interacting fermions in one dimension 
using the conventional linear LL theory. This theory 
is universal in that it predicts the low-energy proper- 
ties of microscopically very distinct systems using only 
very few measurable parameters^ One of its corner- 
stones is the realization that the elementary excitations 
of one-dimensional interacting fermion systems are collec- 
tive bosonic charge and spin density waves^^ This is in 
stark contrast to higher dimensional systems, where the 
Fermi liquid theory^ predicts that the elementary excita- 
tions are only weakly affected by interactions and remain 
fcrmionic. Spin and charge density waves can be encoded 
into the bosonic fields (f> v (x) (y = ,s, c) and the canoni- 
cally conjugate fields 6 v (x), wherei 

[<p v {x),d v 9 ll {y)\ = i-KS^ u S(x - y). (4) 

In the low-energy regime, the spectrum e(fc) of the phys- 
ical fermions can be linearized around the two Fermi 
points, £i?.L(fc) ~ vp(±k — kp). For this spectrum, the 
kinetic energy becomes quadratic in </>„ and Q v . More- 
over, the interaction energy is generally quadratic in the 
charge and spin densities p c . s , 

Pc,s(x) = V<f) c Jx). (5) 

Therefore, the interacting Hamiltonian, despite being 
quartic in fcrmionic operators, remains quadratic in the 
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A(k,u> ^ 0) 


(2n - < fc < (2n + l)k F 




-5( 


' (2n + l)V^c ! <5f 
x/2 


2n ) 


* if i tf-^y 2 


A(k,u> « 0) 


k « (2n± l)fc F 








S(fc,w) 


2nfeF < fc < 2(n + 1)/cf 


u DSF 


1 1 

2 ~ 2 


I V2 2n ) 


2 i 





TABLE I: Exponents for the spectral function A(k,uj) (see Fig.[7]for notations) and the charge density structure factor S(k,uj) 
(see Fig.[8]for notations) at the edge of support. The exponents are determined in terms of the phase shifts 5± = AS± c (k~2nkF) 
and 8± = AS± c [(2n + — fe] calculated in Eq. 1)48 [I. Moreover, m± = (n + 1/2 ± 1/2) mod 2. The exponents and the edge 
of support for the spin structure factors S (k,u)) and S zz (k,uj) coincide with the ones for S(k,uj). 



bosonic basis. It can be shown that the full Hamiltonian 
H becomes a sum of commuting harmonic charge and 
spin terms, Ho = H c + H s , which are given by£ (using 
h=l) 



I 1 



(6) 



for v = c, s. Interactions between the fermions lead to 
different velocities of charge and spin modes, v c ^= v s , 
and thus remove the degeneracy which is present in the 
nonintcracting system. This leads to the spin-charge sep- 
aration that has been observed in experiments^— once a 
spinful fermion is injected into the system, it creates spin 
and charge density waves propagating at different veloc- 
ities, v s and v c . Moreover, interactions determine the 
value of the constants K c s which characterize the dy- 
namical correlation functions and thermodynamic prop- 
erties like the compressibility and the magnetic suscep- 
tibility. We shall assume that the interactions are repul- 
sive. In the charge sector, this entails < K c < 1, where 
K c = 1 corresponds to noninteracting fermions. Further- 
more, repulsive interactions lead to v s < v c . In addition 
to Hq, the spin part of the Hamiltonian will generally 
contain a sine-Gordon term 
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(27m) 2 



dxcos[2v^0s(x)], 



(7) 



where a is a short-distance cutoff. In terms of the 
original fermions, H g corresponds to spin-flip scatter- 
ing cx V , L^'0-£4'0i{^'0-RT an d 9 is proportional to the 2kp- 
componcnt of the interaction potential. A renormaliza- 
tion group (RG) analysis shows that for repulsive inter- 
actions (K c < 1) the interaction strength g flows to zero 
as the bandwidth is reduced, so H g is irrelevant^ In the 
absence of a magnetic field, the system Hamiltonian must 
conserve SU (2)-symmetry, so the components S ^ and 
S zz of the spin density structure factor in Eq. ([T]) must 
coincide. For g = 0, this requirement leads toZ^ 



A"., 



1. 



(8) 



For energies small compared to kp /2m, where m is 
the bare mass of the fermions, all interaction processes 
will involve particles close to the two Fermi points ±kp. 
The physical fermion operators ip a (f =tj i) can then 



be projected onto linearized bands of right-moving and 
left-moving fermions ty aa (x) (a = R, L) with momentum 
close to ±k F , i.e., ip a (x) = e lkpX ipR a (x) + e~ lkFX if; L(J (x). 
These are related to the bosonic fields 9 V and <j) v via the 
bosonization idcntity^iS 



ip ar7 (x) cx 



1 



-i/V2[a<j> c (x)-6 c {x)+acr<p s (x)-<79 s (x)] 



(9) 



The cutoff length a is used to regularize the ultravio- 
let behavior of the theory. One can choose a = 1/A, 
where A <C kp is the width of the linearized bands. In 
principle, the bosonization identity also contains Klein 
factors to ensure the correct fermionic anticommutation 
relations, e.g., {ip aa (x) 7 ip\ T {y)} = 5 a p5 <TT 5{x - y). How- 
ever, we did not write them out explicitly because they 
commute with 6 V and <p v and drop out whenever Eq. @ is 
used to calculate expectation values of an operator which 
conserves charge and spin. The spectral function ([2]) 
contains the expectation value t), - ct(T (0, 0)}) 

which obviously satisfies this requirement and the same 
is true for the other dynamic response functions we shall 
calculate. 

The spectral function of a linear LL can be obtained 
by expressing the chiral fermions ip aiT in terms of <j> v 
and V using Eq. © and then calculating the bosonic 
expectation values with respect to the Hamiltonian ([5]). 
However, the bosonic language is not well suited to in- 
clude band-curvature effects which generate anharmonic 
terms (in boson creation-annihilation operators) in the 
Hamiltonian. Bosonic perturbation theory in these an- 
harmonic terms leads to divergences^ The physical rea- 
son is that, e.g., the right-moving bosonic spin/charge 
excitations with momentum k always have an energy v v k 
[y = s, c), which remains linear in k, so the boson velocity 
v v for right-movers is momentum-independent. Roughly 
speaking, bosons with different momenta thus have a long 
time to interact and this leads to a breakdown of pertur- 
bation theory. It turns out to be beneficial to rephrase 
the problem in a fermionic language by introducing lcft- 
and right-moving fermionic spin and charge quasiparti- 
cles, spinous and holons. The fermionic spectra will be 
curved, thereby avoiding the problem of the bosonic the- 
ory. 

For spinless fermions, the transformation between 
the physical, interacting fermions and noninteracting 
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fcrmionic quasiparticles can be performed directly using 
a unitary transformation^ or via bosonization and sub- 
sequent refcrmionization.— In the spinful case, we use the 
latter option. First, Hq is diagonalized by the canonical 
scaling transformation 



(10) 



Each of the Hamiltonians H v [6 v ,<j> v ] (y = c, s) has the 
form of a free bosonic Hamiltonian with Luttinger param- 
eters equal to unity and should therefore be representable 
as a noninteracting fcrmionic Hamiltonian by reversing 
the bosonization procedure. Indeed, the bosonization 
identity ((9]) can be used to refermionize the theory by 
defining new fermionic operators as 




ip av oc 



1 



-i(acj> v -e v ) 



(11) 



for a = R, L = +, — and v = s, c. Written in terms of 
spinons tp as and holons ip ac , the spin and charge parts of 
the LL Hamiltonian ([5]) assume the form of free fermionic 
Hamiltonians with linear spectrum, 



H v = -iv v ^2 a J dx : ffi av (xjS7i) av (x) 



(12) 



a=R,L 



The colons denote normal ordering with respect to the 
ground state, which is given by filled Fermi seas: for 
a = R (a = L) all states within the bandwidth A 
with negative (positive) momentum arc singly occupied 
whereas all other states are empty. 

The bosonic fields </>„ and 9 U are related to the densities 
of the fcrmionic quasiparticles. In analogy to the conven- 
tional bosonization result^ the refermionization formula 
(fTTT) leads to 



Pav{x) =■ 4>i„(x)ijj a v(x) 



-— V[a&,(x) -9 u (x)}. 
Zn 



(13) 



This allows us to establish the relationship between the 
physical fermions and the spinons and holons. For this 
purpose, the physical fermions tp aa are expressed in terms 
of the rescaled bosonic fields. Using © and (flU)) it is 
straightforward to show that 



i) a \{x) oc ■ip ac {x)F ac (x)ip as {x)F as (x), 
ip a l{x) oc 4> ac {x)F ac (x)ipl ls (x)Fl s (x). 



(14) 



The Klein factors were discarded for the same reason as 
before. The unitary string operators F av (x) are functions 
of the quasiparticle densities and are given by 



Fau(x) = exp < — ia 



dy [8 (y) + S (»)] 

1 

(15) 



where the effects of the interactions are contained in the 
-Kc-dcpcndcnt phase shifts 

S+c 

2ir 

S-c 

2tt 

S+s 

2tt 




= 0. 



(16) 



For the phases S± s , we used K s = 1 which is a con- 
sequence of iS'J7(2)-symmetry As an example consider 
the creation of a right- moving spin- up fermion ifi^- Ac- 
cording to Eq. (fl"4"|) . it corresponds to the creation of a 
holon ip Rc and a spinon ij) Rs as well as of low-energy spin 



and charge density excitations F Rc , F Rs . Conversely, the 
creation of a spin-down particle corresponds to the anni- 
hilation of a spinon. 

For a linear spectrum, the dynamics of the string op- 
erators (|15| is governed by the noninteracting Hamilto- 
nian Hq and becomes very simple. Hence, A(k,u>) can 
be calculated from the refermionized operators and one 
recovers the well-known LL result which is applicable for 
k w zfcfcfi 12 i 13 ' 52 Note that in the conventional approach 
to calculating A{k,oS), one first evaluates the correlators 
as functions of x and t, and only then performs a Fourier 
transform. The latter step is a quite complicated problem 
in contour integration due to the existence of singulari- 
ties in the integrand at x = ±v v t. Our method yields the 
threshold exponents in a much simpler way and provides 
a clear interpretation of A(k,u>) in terms of holons and 
spinons. 

As an example, let us discuss the spectral function for 
u) < and k w +kp- Expressing ip R ^(x") using Eqs. (fT4]) 
and (|15p and Fourier transforming allow us to rewrite 
A(k, uj) as a convolution of correlation functions involving 
right- and left-moving spinons and holons. As k rj +kp, 
we measure the momentum relative to the right Fermi 
point, and use q = k — kp. We then have 

A{q,u) = dxdte**-*-* ^ Rt (x,t)^ m (0,0)) 

1 f dwRc dqFtc duj Lc dq Lc du> Rs dq R , 
~ 2tt 



2tt 2tt 2tt 2tt 

x G Rc (q Rc ,u Rc )GL C (qL C ,ujL C ) 

x Gr s { ( 1Rs^Rs)Gls{ ( 1Ls,^Ls), 



2tt 2tt 



(17) 



where q Ls = q— qRc~ qLc~ <lRs and ui Ls = uj-uj Rc -uj Lc - 
uirs due to energy and momentum conservation. We used 
the Fourier transforms of the correlation functions 



G Rv {x,t) 



G Lv (x,t) 



< e - t 5+„ f» dy~p R „{y,Q)^ Rv ^ (18) 
e iS -" / X d VPLu(v,t) e -iS-v 1° dyp Lv (y.0)\ 



For linear spectrum, the time-dependence of the chi- 
ral left- and right-mover densities is simple, p av {x,t) = 
Pav{x — otv v t). These correlation functions can be cal- 
culated most easily in the bosonic basis, but a direct 
calculation in the fcrmionic basis is also feasible^ For 
the Fourier transforms, one finds 

("+c -A 2 1 

G Rc (q,uj) oc 5(cj - v c q)9(u) + v c q)(u> + v c q)\ 2 " / 

GLc(q, oj) oc 6(lo + v c q)6(u} - v c q)(uj - v c q)( 2w ) \ 

G Rs (q,uj) oc 5(u) - v s q)9(u> + v 8 q)(u} + v s q)\ 2 " / 
G Ls (q, w) oc 8(w + v s q)S(iLi - v s q). (19) 

The arguments of the <5-functions reflect the linear spec- 
trum of the quasiparticles. The distinct form of G Fs {q, w) 
is a consequence of <5_ s = 0. The Heaviside-6* functions 
appearing above have a simple interpretation in terms of 
the Fermi seas of spinons and holons. For right-movers 
(left-movers) all states with negative (positive) momenta 
are filled, thus only excitations with positive (negative) 
momenta can be created. The exponents occurring in 
these correlation functions can be interpreted in the con- 
text of the theoretical treatment of the Fermi edge sin- 
gularity problem by Schotte and Schotte^ 

The quadratic "Anderson" -terms 6+ c , S^ c , and 5+ s in 
the exponents of Eq. ([HI)) are a consequence of the or- 
thogonality catastrophe^ They reflect the shake-up of 
the right-moving and left-moving holons and of the right- 
moving spinons near the Fermi points, respectively, as a 
reaction to the introduction of the hole into the system. 
Left-moving spinons are not shaken up because SU(2)- 
symmetry leads to <5„ s = 0. The exponents linear in <5 +c 
and 5+s, on the other hand, are "Mahan" terms^£ which 
arise because the injection of a right-moving hole leads to 
the creation of a right-moving holon and a right-moving 
spinon. 

Consider A{q,uS) for u> ~ v c q. The singularity at this 
energy is generated by points in the integrand where 
uj w ujrc = v c qfj c . Due to momentum and energy conser- 
vation, the energies wi c and ujr s as well as the momenta 
qLc and qn s will be small, so the functions Gl c and Gr s 
in the integrand lead to power-law singularities. One 
finds the exponent 



/'- 



-1 



2 71 



S+s 

2ir 



(20) 



Similar arguments allow the calculation of the exponents 
at the spinon mass shell. Because v s < v c , the spectral 
function vanishes for \u>\ < v s \q\. It is characterized by 
power-law singularities at the mass shells of spinons and 
holons, A(q,uj) oc (lj— Vc^q)^^- , with the exponent o 12 ' 13 



li- 

/<■: 



1 

1 



- 2 



(21) 




v c q 



v s q 



A{q,uS) 




-v s q -v c q 



FIG. 2: (Color online) The spectral function A(q,uj) for a 
linear Luttinger liquid with repulsive interactions along a cut 
for fixed q — k — kF < 0. A(q,uj) is characterized by sharp 
power-law singularities at the holon and spinon mass shells 
w = v c . s q and at the inverted holon mass shell uj = —v c q. 



In addition, one finds a singularity at the inverted holon 
mass shell, A(q,u) oc (u + Vcq)^^*, where 



1 



Kr 



1 



- 2 



(22) 



The corresponding power-law at ui = —v s q is suppressed 
because the particular form of GLs{qLs,ULs) (|19p entails 
qL S = u!l s = 0. This is a consequence of K s = 1. It 
means physically that the right-moving hole with energy 
oj and momentum k cannot excite left-moving spinons. 
Figure [2] shows a cross-section of the spectral function 
A(q, ui) along a fixed q = k — kp < 0. 



A. Corrections to the linear LL theory for 
nonlinear fermionic spectrum 

If the fermionic spectrum is not linear, two types of 
corrections have to be added to the Hamiltonian Hq (|12|) . 
First, the spectrum of the holons and spinons e c , s {k) will 
start to deviate from its linear behavior. In the case of the 
holons, the leading correction for k — > kp is quadratic, 



1 



e c (k) = v c (k - k F ) + - — -(k - k F ) 
2m* 



(23) 



where m* is an effective mass which is related to the 
bare mass m and the compressibility of the system and 
is generally positive, see Eq. (|77j) . The spinon spectrum 
e s (k) also becomes nonlinear but its form is restricted by 
S'C/(2)-symmetry. The spin- up/spin-down symmetry of 
the physical fermions translates into particle-hole sym- 
metry of the spinons. Therefore, the leading curvature is 
cubic, 



,(fc) = v s (k - k F ) - C{k - k F ) 3 



(24) 



where in general £ > 0. In addition to causing a band 
curvature of the spinon and holon spectra, a nonlinearity 
in the spectrum of physical fermions also leads to interac- 
tions between spinons and holons as well as within these 
two species. The different forms of holon and spinon mass 
shells |23l) - (IMl) entail strikingly different consequences as 
far as the relevance of these interactions is concerned. 
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For the calculation of the spectral function A(k,uj) 
near ui e c (fc), one needs to consider the scattering phase 
between a holon at momentum k and holons and spinons 
at the Fermi edge. To lowest order in the interaction^ 
the scattering phase with left- and right-moving holons 
and spinons is proportional to V° u (k — fcj?)/(wd — av„), 
where V^ u (k—kp) for a = R, L = +, — and v = c, s is the 
corresponding interaction potential and Vd = de c (k)/dk 
is the velocity of the holon with momentum k. For 
k — > kp, the Hamiltonian of the system should reduce to 
that of a linear LL. In particular, this means that spinons 
and holons become noninteracting, i.e., V£ V (Q) = 0. Due 
to symmetry the expansion of the interaction potential 
starts with the quadratic term, V£ v (k — kp) oc (ft — fc^) 2 
for k — > fci?- For a quadratic spectrum, Vd — v c oc (k — Hf) 
for k — > kp, whereas V& ± v s and Vd + v c remain finite 
in this limit. Therefore, all scattering phases vanish and 
interactions of the holon with other quasiparticles do not 
modify the spectral function for k —¥ kp. 

The quadratic form of the holon spectrum (|23p . on 
the other hand, does lead to a change in A{k,u). Let 
us focus again on k < kp and ui < e c (fc) < 0. As de- 
picted in Fig. [31 the injection of a hole with momen- 
tum k and energy to leads to the formation of a holon 
with shifted momentum k + Ak on mass shell and low- 
energy spinons and holons near the Fermi points with 
total momentum Ak. If uj is close to the mass shell, i.e., 
|w-e c (ft)| < \e c (k)-v c (k-k F )\ = (ft- fc F ) 2 /(2m*), en- 
ergy and momentum conservation enforce Ak <C (kp — k). 
Therefore, the "deep" holon at momentum k + Ak be- 
comes well separated from excitations near the Fermi 
points and the respective regions of the spectrum can- 
not overlap, see Fig. [3] In this case, the averages over 
the fermion operator tpfi c (which creates the deep holon 
at momentum k) and the string operators (which create 
the low-energy excitations near the Fermi points) in the 
definition ([15} of Gn c (x,t) can be separated, similar to 
the case of spinless fermions^ 2 - Then, one finds a modi- 
fied exponent near the holon mass shell. For uj m e c (fc), 
A(k, uj) oc [ui - e c (ft)] _M °.- , where 




(25) 



This exponent can again be interpreted in the language 
of the Fermi edge problem: the quadratic exponents S± c 
and 8%_ s are Anderson terms indicating a shake-up of the 
right- and left-moving holons as well as the right-moving 
spinons near the Fermi points, respectively The linear 
exponent S +s can be interpreted as a Mahan term due 
to the creation of a spinon near the right Fermi point. 
In contrast to Eq. ([2TJ)) . there is no Mahan term associ- 
ated with 8 +c because, as explained above, the holon at 
momentum k + Ak is in a different part of the spectrum 
than the Fermi point. 

This exponent differs from the LL result ([2TH) . The ex- 



v,.Ak 



\e D (k + Ak) - ui\ 



Ak 



Ak 



tc(k)/ / 



k 



k f 



(k,uj) 



k F ) 



FIG. 3: (Color online) The injection of a hole with momentum 
k and energy uj ~ e c (k) leads to the formation of a holon on 
the mass shell and low-energy excitations around the holon 
Fermi edge. For clarity, the created low-energy spinons are 
not displayed. For this particular configuration, energy and 
momentum conservation lead to v c Ak — \u — e c (k + Ak)\. The 
impurity band (around k) and the low-energy band (around 
k F ) can be separated only if Ak <C k F — k. This is the case 
for \uj - e c (fc)| < (k - fc F ) 2 /(2m*). 



ponent _ holds in a region of width {k — kp) 2 / (2m*) 
around e c (k). Beyond this window, the curvature of the 
spectrum becomes irrelevant and the exponent crosses 
over to the LL exponent /il . Note that even for K c — > 1 , 
the exponent /Ltg _ — > \/2 — 1 differs from the LL predic- 
tion 1/2. 

For the calculation of the exponent near the spinon 
mass shell, uj f=a e s (fc), the interactions among spinons 
become important. This is a consequence of the cubic 
term in the spinon spectrum ([24| . Consider the scat- 
tering of a spinon with energy uj on low-energy left- 
and right-moving spinons and holons near the Fermi 
points. As previously, the lowest-order expansion of the 
respective interaction potentials for k — > kp is quadratic, 
V a v{k — kp) oc (k — kp) 2 . The velocity of the spinon is 
given by Vd = de s (k)/dk. However, in contrast to the 
quadratic holon spectrum, the leading curvature of the 
spinon spectrum e s (k) is cubic, so Vd — v s oc (k~kp) 2 . As 
a consequence, the scattering phase V+ S (fc — kp) / (vd~v s ) 
remains finite even in the limit k — > kp. Hence, the 
scattering among spinons cannot be treated as a small 
perturbation. We shall investigate the scattering phase 
shifts for excitations near the spinon mass shell in the 
next section. 



IV. QUANTUM IMPURITY HAMILTONIAN 

The investigation of the spectral function for mo- 
menta away from ±kp necessitates a comment about 
the momentum conservation when decomposing an in- 
jected physical particle or hole with momentum k into 
a spinon-holon pair. Let us start from noninteracting 
fcrmions with Fermi momentum kp. In terms of mcasur- 
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FIG. 4: Holon spectrum uj c (k c ) and spinon spectrum uj a (k s ) 
in the noninteracting case. The holon spectrum has Fermi 
momentum 2kp. The spinon spectrum is defined for — kp < 
k s < kp and the Fermi point is placed at k a — 0. The spinon 
spectrum is particle-hole symmetric, so the dashed line indi- 
cates energies of spinon holes. 

able quantities, kp can be defined as the smallest positive 
momentum k for which A(k,uj = 0) 7^ 0. According to 
this definition, kp is defined by the singularities in the re- 
tarded Green's function, which according to Luttinger's 
theore m 56 ' 57 are not shifted when the interactions are 
turned on. Therefore, the value of kp is not affected by 
interactions. 

The charge density of the physical fcrmions is related 
to kp by p c — 2kp/ir, where the factor 2 results from the 
two spin orientations. The spinful physical fcrmions can 
be expressed in terms of two species, spinons and holons, 
of spinless fermions. Since only the holons carry charge, 
their density must be equal to the physical charge den- 
sity. Hence, the holon Fermi momentum k F = 2k p. This 
has been found in the case of a generic strongly inter- 
acting system 5 ^ as well as for intcgrable models at any 
interaction strength ! 59 ' 60 In particular, the Bethe ansatz 
solutions for the Hubbard model at low filling in both the 
noninteracting limit and in the limit of infinite interac- 
tion lead to a parabolic holon spectrum ui c (k c ) shown in 
Fig.H 

In the same context, it has also been shown that 
the spinon momentum is only defined up to multiples 
of 2kp. This is well illustrated by the limit of strong 
repulsion^ In this limit, spinons live on a lattice pro- 
vided by the holons. Therefore, we restrict the spinon 
spectrum uj s (k s ) to momenta — kp < k s < kp and place 
the spinon Fermi level at zero momentum, w s (0) = 0. 
For a physical particle with momentum k, the momen- 
tum conservation can then be expressed as^ 

k = k c + k s ±k F for k ^ 0. (26) 

The spinon and holon mass shells e c (fc) and e s (k) can be 
constructed from the spectra u) c (k c ) and uj s (k s ), respec- 
tively, by a combination of shifts and inversions. 

The spectral function of interacting one-dimensional 
systems is generally characterized by power-law singu- 
larities. To explain the physical mechanism, let us focus 
briefly on the spinless case and consider a spinless hole 
with momentum k < kp tunneling into the system. If 
the energy is on the mass shell, the final state will con- 
tain a single "deep" hole with momentum around k and 
energy e(fc). If the energy is close to, but not precisely 
at, the mass shell, the final state will contain the deep 



hole and, in addition, a number of particle-hole pairs, 
see Fig. [3] Since the energy available for the formation of 
these pairs is small, they will be located close to the Fermi 
points ±kp . It has been shown perturbativel y 20 ' 24 that 
for a nonlinear excitation spectrum the deep hole can be 
regarded as separate from the particle-hole pairs at the 
Fermi points. The Hamiltonian can be projected onto 
three subbands, one containing the deep hole and two 
others containing excitations near the two Fermi points. 
The functional form of the spectral function is deter- 
mined by how the deep hole interacts with the particle- 
hole pairs at the Fermi points. The physical mechanism 
is thus similar to the Fermi edge problem ; 47 ' 48 ' 55 albeit 
in the present case with a mobile impurity instead of a 
static scattering potential. 

The previous argument extends to spinful systems: in 
this case, power-law singularities occur at the spinon 
mass shell e s (k). Let us consider the case of a right- 
moving spinful hole with momentum < k < kp and 
energy close to e s (k). The final states giving rise to the 
power-law singularity will contain a deep spinon impurity 
with momentum k s = k — kp =: k c i < 0, a holon near 
its Fermi momentum and additional spin-carrying and 
charge-carrying particle-hole pairs at all Fermi points. 
Once more, the system can be projected onto narrow sub- 
bands containing the impurity and the Fermi points, re- 
spectively. This corresponds to decomposing the spinon 
operator by retaining only Fourier components close to 
zero and to kd, ipR S — > ipRs + e x d. The spectrum of the 
states near the Fermi points can be linearized and they 
are thus described by the LL Hamiltonian, 

Hp = -i v v a J dx : ^(a^VVw^) : . 

v=c,s a=R,L 

(27) 

Located on the spinon mass shell, the impurity has a 
velocity Vd = de s (k)/dk. The Hamiltonians containing 
its kinetic energy and the interaction between the mobile 
impurity and the subbands at the Fermi edges are given 
by 

H d = J dx dt(x)[e s (k) - iv d V]d s {x), (28) 
H mt = J dx^2 V au (k)p av (x)dl(x)d s (x), (29) 

where p au (x) =:^ av {x)ip av (x) : denotes the quasiparticle 
density in the bands around the Fermi edges, and the 
interaction constants V au {k) are yet to be determined. 

In addition, the interacting system Hamiltonian will 
generally contain spin-flip interaction terms of the form 

Y, 4(P)^(P + l)V(q)4(P)MP - l)- (30) 
p,p',q 

The projection of this term onto a reduced band struc- 
ture with bands around the Fermi edges and around the 
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the rotated basis reads 



FIG. 5: (Color online) Spinon interaction process generated 
by the projection of the spin flip interaction in Eq. (|30[) : a 
spinon particle-hole pair near the left Fermi point scatters 
into a spinon hole at k a < and a spinon near — k a > 0. 
The corresponding interaction operator is similar to the sine- 
Gordon term ([7]) and, analogously, its amplitude flows to zero 
upon bandwidth reduction. 



impurity state leads to the following terms: (i) a density- 
density type interaction between particles near the Fermi 
edges X Such a contribution merely renormalizes the Lut- 
tinger parameter K c , as well as the spinon and holon ve- 
locities, v c and v s ; (ii) a density-density type interaction 
between the impurity and particles near either the left 
or right Fermi point. This leads to a term of the form 
(|29|) and can be absorbed into V a „; (iii) spin- flip inter- 
action between particles at the two Fermi edges. The 
bosonization of these terms leads to a sine-Gordon term 
([7]) which vanishes logarithmically when reducing the 
width of the linearized bands around the Fermi points. 
In the SU(2)-symmetric case, K s approaches unity simul- 
taneously with the reduction of g in Eq. [7] to zero; (iv) 
spin-flip terms of the type e~^ % ^" +e ^b\d s which trans- 
fer two left-moving spinons with k « to the d s -impurity 
band (k s < 0) and to its mirror image (k — kd, the "& s - 
band" ) . The process is illustrated in Fig. [5] These terms 
are the finite- k counterparts of the sine- Gordon term in 
Eq. (JT)). Similar to the latter, these terms vanish in the 
limit of low energy of excitations. 

Therefore, considering an interaction term of the form 
([2T)| is sufficient. The interacting Hamiltonian (j2"T)) - (|2l)]) 
can be diagonalized using a unitary transformation, 

U = cxp < -i / rf^y^ AS av N av (x)dl(x)d s (x) > , 

{ "'- 00 av J 

(31) 

where N av = J x dyp au {y). The unitary operator U is 
characterized by the phase shifts 



AS av 



vd - av v 



(32) 



This transformation removes the interaction term, 
W(H + H d + H mt )U = H + H d . The d s -operator in 



UU s (x)U = d s (x) 



exp < —i 



dy^2,t±5 av p av (y) 



(33) 



One can see from Eq. (fl"5"j) that A(*> Q „ simply adds to 
the phase 5 au which resulted from the refermionization, 
to produce the total shift b* av = 5 av + A5 au . Power- 
law singularities at the spinon mass shell in all dynamic 
response functions will be characterized by the phases 
<5*„. The hole sector (cj < 0) of the spectral function at 
momentum k is given by 

A(k, u<0) = — [ dtdxe tkx - lut 
2tt J 

H„+H d +H int , (34) 

and is independent of a. It can be calculated 
by representing the physical fermions as i^Rf ~ 
e lkx ip Rc Fn c d s Fn s . The correlation function of the impu- 
rity after application of the unitary transformation (|31[) 
is 

(dl(x,t)d s (0,Q)) Hd =e^ k)t 5(x-v d t), (35) 

and the correlation functions involving the string opera- 
tors can be calculated by using the bosonized expressions. 
One finds that A(k,u) cx [uj — e s (fc)]~ At o>- for uj w e s (k), 
with the exponent 



2tt 



2n 



(36) 



We have seen previously that the cubic spectrum of the 
spinons means that the phase shift A<5 +s remains finite 
even for k — > kp and cannot be treated as a perturba- 
tion. Luckily however, S'J7(2)-symmetry can by used to 
fix 6± s = 5± s + A(5± s for arbitrary momenta. As we shall 
show in Sec. IV G[ the exponents of the power-law sin- 
gularities in the spin correlation functions S zz (k,ui) and 
S ^(k,cj) can be calculated as functions of 5± s . In an 
SU (2)-invariant system both exponents have to coincide 
for all momenta and this is only the case for 



J ±s 



0. 



(37) 



As we laid out after Eq. (|T5|) . the phase shifts 5± 8 can 
be interpreted as the probability amplitude of a shake-up 
of the spinons near the Fermi points. Vanishing phases 
5± s = mean that the creation of the spinon d s due to 
the incoming hole does not lead to such a shake-up and 
thus to a formation of particle-hole pairs in the spinon 
sector. For \u\ > |e s (fc)|, any excess energy can only 
be used to create particle-hole pairs in the holon sector. 
In this important sense, spin-charge separation remains 
meaningful even for systems with band curvature. 
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V. DYNAMIC CORRELATION FUNCTIONS 
AWAY FROM THE FERMI POINTS 

A. Phenomenology for the scattering phase shifts 



In the limit k — 
latcd from Eq. (|36 



kp, the edge exponent can be calcu- 
by using AS± C — in addition to the 
relation 8*L S = which is valid at arbitrary k. It turns out 



that /Xq_, unlike 
LL exponent 



coincides with the corresponding 



Mo, 



(38) 



Therefore, for u> w e s (k), one finds A(k,ui) oc [lu — 
e s (/c)] _M - even in the presence of band curvature. 

In this section, we shall extend this result to momenta 
k away from the Fermi points. We argued that for repul- 
sive interactions, the edge of support of A(k as +kp,u>) 
coincides with the spinon mass shell u> « e s (fc), and at 
this edge A{k,w) is characterized by a power-law singu- 
larity. The state responsible for this singularity contains 
a spinon on mass shell and a holon at the Fermi point, 
so the mobile impurity d s has the quantum numbers of 
a spinon. Since the edge exponents must be continuous 
as a function of k, this must still be true for momenta 
away from the Fermi points. Hence, in order to calcu- 
late the dynamic response functions in this region, we 
can still make use of the same mobile-impurity Hamil- 
tonian but in contrast to the limit k — » kp, we can no 
longer use the lowest-order expansion of the phase shifts 
AS± C . Instead, we shall derive phenomenological expres- 
sions which relate the phase shifts to measurable prop- 
erties of the spinon spectrum e s (k). In addition to the 
possibility of directly measuring e s (k), this function can 
be evaluated numerically by well-developed routines. It 
can also be calculated exactly for integrable models. 

Let us focus again on the case k < kp and investi- 
gate the spectral function in the vicinity of the spinon 
spectrum uj ss e s (fc). The configuration responsible for 
the edge singularity contains a spinon impurity near 
kd = k — kp, a holon close to the Fermi point, as 
well as particle-hole pairs of spinons and holons at the 
Fermi edges. The effective Hamiltonian is given by 
Eqs. (j2"T)) - (f2T)|) . We express the quasiparticle densities 
Pau{x) in terms of bosonic operators Q v and </>„. The 
interaction between the impurity d s and the low-energy 
spinons p as vanishes because £j_ s = 0. The Hamiltonian 
H = Hq + Hd + Hi n t can be written as 



ffo = 7^ / dx 
Zn 



K C (V6 C ) 2 + — (V0 C ) 2 



^-Jdx [(V0 S ) 2 + (V0 S ) 2 ] , 
Hd= I dx d\(x)[e s {k) — ivdV]d s (x), 



H^t = I dx 



V R V- 



2- 



- VlV- 



2tt 



dUs, 



(39) 



and contains two parameters, Vb(fc) and VrQc), which de- 
scribe the interaction between the impurity and holons 
at the Fermi points. The interaction term can be re- 
moved using a unitary transformation like Eq. pip and 
leads to the following relations between phase shifts and 
interaction constants, 



^Kc(V L + V R ) 
Vl-Vr 



-A6- C (v d + v c ) - AS +c (v d - v c ), 

-AS- c (vd + v c ) + AS +c (v d - v c ). 

(40) 



In order to fix Vl and Vr, we need two relations. The 
first one can be derived by considering the response of 
the system to a uniform charge density variation. Since 
the spinon spectrum e s (k) is defined with respect to the 
chemical potential [i, for fixed k the shift in the edge 
position upon a variation of the density 8p c is given by 



5E 



de s (k) 
dp c 



dp c 



Sp c = 



de s (k) 
dp c 



IK, 



8p c 



(41) 



In the last equality, we used a phenomenological relation^ 
between K c and the compressibility k = dp c /dp = 
2K C / (irv c ). A second way to calculate the same energy 
shift is to use the Hamiltonian (|3U)) . According to Eq. © , 
a density variation leads to a finite expectation value 
(V^c) = —TrSpc/y/2. We use the Hamiltonian (f3"9"f to 
calculate the energy of a state containing a spinon impu- 
rity d s with momentum kd and a holon at the Fermi edge. 
Then, we investigate the change of this energy due to the 
density variation. In this case, the interaction term (|39[) 
shifts by 



Vr + Vl /V7 ^ . Vr + Vl , 

OUint = " (V0c) = 7= 0p c . 



2n 



2V2 



(42) 



The impurity momentum is given by kd = k — kp. As 
the density variation affects the Fermi momentum, Skp = 
ir5p c /2, one obtains an energy shift in the operator Hd 
by 



SE d = - 



nvdSpc _ 7T de s (k) 
2 _ _ 2 dk 



5p c . 



(43) 



Finally, the energy of the holon at the right Fermi point 
changes due to the shift of the chemical potential and 
leads to SEq = (dp/dp c )8p c . Comparing the shift (|4Tj) to 
the shifts 5Eo + SE d + SEi n t calculated using the effective 
Hamiltonian, one obtains the first relation, 



V R + V L de s (k) 



2V2 



dp c 



7T de s (k) 
2 dk ' 



(44) 



In order to obtain a second relation which fixes the differ- 
ence Vl — Vr, we consider the effect of uniform Galilean 
boost of the system. Let us assume the liquid is moving 
at velocity Su and the incoming particle has momentum 



12 



k. In a reference frame moving with the liquid, the mo- 
mentum of the injected particle is k — mSu, where m is 
the bare mass of the physical fcrmions. The total change 
of the energy acquired by the liquid due to the boost is 
given by 



5E' = — 
2m 

= Su 



(k — mSu) 2 



2m 



,{k — mSu) 



de s (k) 



dk 



(45) 



Next, we calculate the same shift using the Hamiltonian 
(f3T))) . The effect of a boost with velocity 5u on the phys- 
ical fermions is to shift the Fermi momentum of right- 
and left- movers, kp' L — > kp ± mSu. This gives rise to 
a difference between right-mover and left-mover density, 
PR ~ Pl — 2m6u/iT. This density difference corresponds 
to a finite expectation value (V0 C ) = y2mSu. We use 
again the Hamiltonian (f3"9")l to calculate to energy of a 
state containing an impurity d s and a holon at the Fermi 
point. Now, we investigate how this energy changes due 
to a finite (V0 C ) = y/2m6u. In the interaction Hamilto- 
nian Hi n t, one obtains 



SEL 



Vl - Vr 
2tt 



(V0 C ) = V2m6: 



V L -V R 



(46) 



The energy shift due to the holon at the Fermi point is 
given by 6E' = K c v c mdu = vpmSu. Finally, the mo- 
mentum of the impurity kd changes due to the shift of 



kp and leads to SE', 



8E' + 5E> with 



—mSuvd- Equating the shifts 
451) leads to 



Vl-Vr 
2tt 



V2r. 



(47) 




Equations (|44[) and (|47j) allow us to fix the interaction 
strengths Vl and Vr in terms of the derivatives of e s (k) 
with respect to the density and the momentum. These, 
in turn, can be related to the phase shifts using Eq. (|40)) . 
The result is 



(48) 



and it is valid for all systems with Galilean-invariant mi- 
croscopic interactions. These relations are valid for all 
— kp < k < kp. The knowledge of these phase shifts 
allows the calculation of all dynamic response functions 
for energies close to the spinon mass shell at arbitrary 
momenta. Note that for the case of the exactly solvable 
ID Yang-Gaudin model ; 61 ' 62 the phase shifts predicted 
by Eq. (|48)) coincide with the exact results obtained us- 
ing the Bethe ansatz (which can be obtained as a limiting 
case of ID Hubbard model considered in 1631). 



B. Edge exponents of the spectral function 

The phase shifts <5* v allow us to derive the edge expo- 
nents of A{k,uj) for momenta k away from ±kp. In the 
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FIG. 6: (Color online) Mobile-impurity band structure for the 
calculation of A(k,u>) for 3/cf < k < 5kp (i.e., n = 2) and 
uj ~ \e s (k — 4/cf)| • The final state contains a spinon with mo- 
mentum kd,2 = k — 5kp < 0, a holon at its Fermi momentum 
as well as two particle-hole pairs in the holon sector which 
absorb the extra momentum 4/cf- In this figure, we assume 
m_ = 2, so the spinon sector contains two additional particle- 
hole pairs. The Hamiltonian is projected onto narrow bands 
around these momenta. Within each band, the spectrum can 
be linearized. 



following, we shall focus on a general momentum inter- 
val (2n — l)kp < k < (2n + l)kp (with integer n) and 
energies near the spinon mass shell, uj ~ e s (fc — 2nkp) 
and its shadow bands, uj ~ — e s (k — 2nkp), see Fig. [7] 
Similar to Scc. lIVl for uj < injecting a hole into the sys- 
tem will create a spinon impurity with momentum near 
kd.n = k~ (2n + l)kp on its mass shell as well as a holon 
at approximately the Fermi momentum. The remaining 
momentum 2nkp must be given to additional spinon and 
holon particle-hole excitations near the Fermi edges, see 
Fig. [5] In terms of the physical fcrmions, the most gen- 
eral configuration of additional particle-hole pairs reads 
as follows, 

$2„ feF := (4>R t T R ' (i>Rl) nai (^Lt)" Lt 0>Li) nLi ■ (49) 

All n aa are integer numbers. Positive values correspond 
to powers of the annihilation operators ip aa , whereas neg- 
ative numbers denote powers of the creation operators 
i\j\ a . Since spin and charge of the incoming hole are ab- 
sorbed by the deep spinon at momentum fc^n and the 
holon at the Fermi point, the total charge and total spin 
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of these additional excitations must vanish. This means 



n-flt + rim + riL-f + n Li = 0, 
nn\ - n R i + n L \ - n Li = 0. 



(50) 



Prom these equations, one finds — ~ur^ and ul\, = 
—tiri- As an additional requirement, the excess momen- 
tum 2nkp must be accommodated. This leads to the 
constraint + — — = 2n and therefore 
to TiRf + TiRi = n. This leaves one free parameter m_, 
which must satisfy the selection rule 

m_ = n(mod 2). (51) 

The general solution for integer n aa reads 



tiri 



n + m_ 
2 

n — m_ 



(52) 



Therefore, the state (|4U)) can be written as 



1>, 



2nk F = Wr^RX^L^Wli) 



f x-( n+TO _)/2 r , / ,t 



Wl^riT • (53) 



Physically, the first term comes about as a result of spin- 
flip scattering. Indeed, the fermionic representation of 
the sine-Gordon term ([7]) produces exactly this type of 
scattering. The second term comes about due to scatter- 
ing between the right and left Fermi points and absorbs 
the excess momentum 2nkp. The spectral function can 
be calculated by decomposing the fermionic operators ac- 
cording to 



ipRf = e ikx il) Rc F Rc d s F Ra x $2nk 



(54) 



and bosonizing 4>2 n fc F and ipR C using Eqs. (j9|) and ([TTjl . 
respectively. One finds that in the interval (2n — l)ftp < 
k < (2n + l)kp, it has a power-law singularity A(k, uj) oc 
[to + \e s (k — 2nkF)W~^ n -~ with the exponent 



_ 1 f (2n + l) y/K, 
2 I V2 




1 



AS 



+c 



A<S_ C 



2tt 



m_ 



(55) 



where AS± C = A5± c {k — 2nkp) is evaluated for momenta 
on the main spinon branch, — kp < k — 2nkp < kp. 
Since the selection rule requires integer m_, the leading 
exponent is given by m_ = for even n, and m_ = ±1 
for odd n. 

A similar line of reasoning can be applied to calculate 
the edge exponents /i* + for (2n—l)kp < k < (2n + l)kp 
and oj > (see Fig. [J])- The configuration with this 
combination of momentum and energy contains a spinon 
with momentum kd = k — (2n — l)kp > as well as a 
holon near the Fermi point. Similar to the previous case, 
the remaining momentum 2(n — l)kp must be absorbed 
by additional particle-hole excitations of the form (14^1) . 



Stipulating again charge and spin neutrality leads to a 
modified selection rule m + = n + l(mod 2). As a result, 



one obtains A(k,ui) oc [a; — \e s (k — 2nkp) 
/!* i is similar to Eq. (|55|) , 



.+ where 



_ 1 / (2n + l)V K 
2V V2 




If 1 



2 \V2K~c 



AS +C - A6_ c 
2ir 



-+ , 



(56) 



except for the different selection rule for ni + . 

Using the formula for the phase shifts (|48|) . it is thus 
possible to express the exponents ([5"5j) and (|55|l entirely 
in terms of the spinon spectrum e s (k) and the Luttingcr 
parameters. This establishes a relation between two dis- 
tinct sets of observable quantities, exponents and spectra, 
which applies to a large range of systems and which can 
in principle be checked in experiments. 

As a consequence of the symmetry of the spinon spec- 
trum e s (k) = e s (—k), the edge exponents also satisfy the 
k —> — fc symmetry, /i* ± (fc) = juf_ n ± (— k). Moveover, 
since the phase shifts AS± C vanish for k — > kp, one can 
also verify that the edge exponents change continuously 
when crossing a Fermi point, 

< ± [(2n + l)k F - 0] = /4+i lT [(2n + l)k F + 0]. (57) 

The final result^ for the edge exponents of A(k, ui) in 
all regions of the (k, w)-plane is shown in Table Q] The 
positions of the edges and the corresponding notations 
for the exponents are shown in Fig. [7J 



C. Charge and spin density structure factors 

The charge density structure factor is defined as 
S(k,u) = I dxdte iut - lkx {p c {x,t)p c (0,0)) , (58) 



where p c {x) = 'Yl (T '^] r {x)'4>a{x) denotes the charge den- 
sity. In the case of a linear spectrum, spin and charge 
density waves with momentum k have a uniquely de- 
fined energy v s>c k. One of the predictions of the linear 
LL theory is that even for nonzero interactions, these 
are stable excitations. Therefore, one finds SpL{k,uj) = 
2K c \k\5(uj — v c \k\) where the effect of interactions is lim- 
ited to a renormalization of the velocity v c =^ vp and the 
pref actor is fixed by the /-sum rulei 

In the case of finite band curvature, a one-to-one re- 
lation between momentum and energy of particle-hole 
pairs no longer exists. Let us illustrate the consequences 
for the noninteracting case. The fermion spectrum is 
given by e(fc) = (k 2 — kp)/{2m) and it turns out that for 
< k < 2kp, S(k,Lo) is nonvanishing only in the interval 
u;_(fc) < uj < tu + (k), where uj±(k) = vpk ± k 2 /(2m), see 
Fig. [5] The density excitation of lowest energy w_ for 
given momentum k contains a particle at the Fermi level 
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FIG. 7: (Color online) Spectral function A(k,uj) in the (k,uj)- 
plane (a) and along a cross-section for fixed < k < fc_p (b). 
(a) A(k,uS) is nonzero in the shaded areas. For repulsive in- 
teractions, the edge of support is at the spinon mass shell 
u — ±|e s (fc)|. At the edge for ui ^ in the momentum range 
(2n — 1)/cf < k < (2n + l)kF, A(k,uj) has a power-law sin- 
gularity characterized by the exponent p,^±. A sharp power- 
law at the holon mass shell e c (k) (dotted lines) exists only for 
k Ri (2n ± l)fe_F. The corresponding exponents /4i,- are con- 
tinuous at uj — 0. (b) Compared to the Luttinger liquid case, 
the edge exponents is modified. The singularity at the holon 
mass shell e c (k) is smeared out because the band curvature 
leads in general to a finite holon lifetime. 



and a hole at kp — k (see inset of Fig. [5]). The upper 
threshold w+ reflects the highest energy such an excita- 
tion can have. It corresponds to a particle at kp + k and 
a hole at the Fermi level. 

For spinless weakly interacting systems, a power-law 
singularity develops at w = w_(fc) but the structure fac- 
tor still vanishes below it, S(k,uj) = for u> < uj^(k)^ 
The upper threshold, on the contrary, no longer exists 
since density excitations can give away excess energy 
to create other particle-hole excitations. Instead, for 
lu ^ uj + (k), the density structure factor decays as a 
power-law. A schematic picture of the domain of sup- 
port of S(k,uj) in the (fc,w)-plane is shown in Fig. [5] 

For spinful systems, the charge density structure fac- 
tor describes responses in the charge sector. Its analog 
in the spin sector is the spin structure factor. Since the 
calculations of both functions are identical and the in- 
variance of the spin structure under spin rotations allows 
us to fix the phase shift <5j_ s , we shall henceforth focus on 
the latter. 

In the absence of a magnetic field, the system Hamilto- 
nian is 5'C/(2)-invariant. One of the consequences of this 
symmetry is identical power laws in different components 




FIG. 8: (Color online) Density structure factor S(k,uj). For 
noninteracting systems S(k,Lj) is constant for u-(k) < uj < 
to+(k) and vanishes otherwise. For nonzero interactions, a 
power-law singularity with exponent p„ F forms at the lower 
threshold cj-(k). In the inset: Particle-hole excitation with 
momentum k giving rise to the singularity at uj-(k). 



of the spin structure factor. In particular, the functions 
S~ + (k,u) = J dxdte iut - ik * (S-(x,t)S+ (0,0)) , 
S zz (k,uj)= [ dxdte lult - lkx (S z (x,t)S z (0,0)) (59) 



must be identical even in the presence of band curvature. 
In terms of the physical fermions, the spin density is given 
by S(x) = | J^a a' 4'(j(x)f aa iip cr '(x), where f denotes the 
vector of Pauli matrices. 

To be specific, let us first investigate S h (fc,w) in the 
n = band, i.e., for < fc < 2kp. Similar to the dynamic 
structure factor S(k,Lu), the spin structure factor will 
vanish below a threshold energy w_(/c) and have a power- 
law singularity at this threshold. The excitation which is 
responsible for the singularity contains a spinon impurity 
at momentum k as well as a spinon hole with small mo- 
mentum. Therefore, we can project the spinon operator 
on two bands using ipR S (x) — > ipR S (x) + e lkx d s (x). For 



S R (x), one finds 



S+(x) = f* m {x)ip Ri {x) 



(60) 
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-ikx 
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2tt 



d\ exp 1 i | — 
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V2 



,(x) - 4> s (x)] 



For the z-component of the spin density, we use 
Sr(x) = - ip f m (x)tp m (x) - i^ Ri {x)^ Ri {x) 

= 2jj Rt (x)lp m (x) - \pRc{x), 



(61) 



where p Rc = ip R ^ xpn-f denotes the right- moving 

charge density. If the spectrum of the physical fermions is 
linear, the spin-charge separation ensures that pR C does 
not involve spinon operators. For a nonlinear spectrum, 
this strict separation does not apply any more and pr c 
may indeed contain spinon operators. This fact will be- 
come important in the calculation of the charge density 
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structure factor. Since all terms lead to the same series 
of exponents, it is sufficient to only project the first term 
and use 



S z R (x) <x e~ ikx dlF^ Rs F R8 + h.c 
oc e~ ikx d\ 



(62) 



I + ) 



h.c. 



The ensuing calculation of the exponent of S ^ is per- 
formed using the mobile impurity Hamiltonian (|27l) - (|29[) 
and general phase shifts <5*„ = 5 au + AS al , and results in 
the leading edge exponent (J,qq. 

In analogy to the previous section, the marginally irrel- 
evant spin-flip scattering ([7]) may create a final state with 
identical energy and momentum but which contains ad- 
ditional particle-hole pairs and thus produces sublcading 
exponents. To calculate these, the operators S R and S R 
must be multiplied by the operator &2nk F (|53[) for n = 0. 
A state arising from \m\ spin-flip scattering events will 
contribute the exponent (for m € Z) 



2 V y/2 
2 VV2 
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2tt 
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(63) 
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where 6± s is evaluated at momentum k € [0,2kp] and 
AS± C at momentum hp — k. An analogous calculation 
yields the leading and subleading exponents for S zz , 
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(64) 
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It should be emphasized that each of the exponents in 
Eqs. (p351) - ([6"4"|) actually generates an infinite "ladder" of 
exponents differing by an integer power, and Eqs. (|63[1 - 
(|64p give only the most divergent exponent. Indeed, 
when we calculate, e.g., the exponent of S l ~(£;,w), the 
expected behavior of mth term is 

S-+(k,w) oc [u - \e s (k F - k)\]~^R[u - \e s (k F - k)\], 

(65) 

where R is at most logarithmically divergent for uj — 
\e s (kp—k)\ — > 0. A Taylor expansion of this function gen- 
erates an infinite ladder of exponents. S'J7(2)-symmctry 
requires that these full sets of exponents generated by 
different /1q ^ and /i§ z m should coincide. This is a weaker 
requirement compared to coincidence of sets of /x^"^ and 



/1q z to . However, it turns out that both these constraints 
lead to the same requirement 



0. 



(66) 



By construction, this equality holds for arbitrary mo- 
menta. Note that in order to fix both phase shifts &J. S , a 
comparison of only the leading (m = 0) exponents is not 
sufficient. 

In the momentum interval 2nkp < k < 2(n + l)kp, a 
spinon impurity will be created at momentum k — 2nkp 
and the second spinon near zero momentum. As in the 
calculation for A(k, u)), the excess momentum 2nkp is 
accommodated by scattering across the Fermi points. A 
number |m| of spin- flip events will lead to a selection 
rule of the form (fSTj) . If we focus only on the lead- 
ing (m = 0) exponent, we can conclude that the spin 
structure factor near its edge of support is given by 
S~+{k,u) oc S"Qe,u) oc {w- |e s [(2n+l)/c F -A:]|}-^ SF 
where 



SSF 




(67) 



with A5± c = A6± c [(2n+l)kp — k}. Here, the phase shifts 
are related to the spinon spectrum by Eq. (|48|) . 

For the calculation of the charge density structure fac- 
tor S(k,uj), the operator p c {x) needs to be examined. 
For a nonlinear spectrum of the physical fermions, this 
operator will contain spinon operators. Therefore, p c can 
create a state which contains a spinon with momentum 
k on its mass shell and additional spinons and holons 
near the Fermi points. For a given momentum k, this 
state is the one with the least energy, so the edge of sup- 
port of S(k,uj) coincides with the one for 5* and 
S zz (k,ui). Therefore, we project p c onto two subbands 
using again ip Rs (x) — > tp Rs (x) + e lkx d s {x). Because the 
respective terms in the projection of the operator p c (x) 
are identical to the terms in the projection of S R (x) in 
Eq. (|62]) . it follows that for 2nk F < k < 2(n + l)k F , 
near the edge of support S(k, uj) oc {lu — |e s [(2n + l)k F — 
fc]|}-^" SF , and 



,DSF 



,SSF 



(68) 



An overview of the edge exponents of the charge and 
spin density structure factors in the different regions of 
the (fc, a;) -plane is contained in Table 

For small momenta k <C kp, the spin and charge 
density structure factors were investigated recently 
within the Abelian^ and non-Abelian^ bosonization ap- 
proaches. In the former approach, the nonlinear spec- 
trum of fermions produces terms which mix the spin and 
charge fields in the bosonized Hamiltonian. The pertur- 
bation theory in these terms developed in Ref. 1 1911 cor- 
roborates one of our conclusions (and of Ref. |18|): the 
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edge of support of S(k,uj) actually coincides with that 
of S zz (k,u). For repulsive interactions, it is located at 
the spinon mass shell, i.e., uj = v s k for small momenta. 
However, the second-order bosonic self-energy actually 
diverges at this edge. As a consequence, it was not possi- 
ble in Ref. [l|| to investigate the shape of the edge singu- 
larity. Despite this drawback, the Abelian bosonization 
provides an easy access to S(k,u) and S zz (k,ui) away 
from the singularities, for instance to the high-frequency 
tails for uj ^> v c k. These "off-shell" tails were studied in 
Ref. [l!| . The broadening of the peaks in S(k, oj) near the 
spinon and holon mass shells was addressed in Ref. fl8j |. 
Its conclusions regarding the threshold exponents agree 
with the limit k — > of our results. 



VI. HOLON EDGE SINGULARITIES 

So far, we have only discussed those power law singu- 
larities which occur at the edges of support of a dynami- 
cal correlation function. We argued that in the case of re- 
pulsive interaction, the spinon-holon excitation with the 
lowest possible energy for a given momentum contains 
a holon at its Fermi edge, whereas the remaining mo- 
mentum and the whole energy are carried by the spinon. 
Therefore, the exponents at the edges of support are char- 
acterized by the phase shifts A5± c produced by the in- 
teraction of a spinon impurity with low-energy holons. 

For energies above the edge of support it is also possible 
to create spinon-holon excitations which contain a spinon 
at its Fermi point and give the entire energy to a holon. 
For a generic system, such an excitation will not be stable 
since energy and momentum conservation allow the decay 
of a holon through the creation of spinon pairs. This will 
lead to a broadening of the threshold. One important 
exception is the case of integrable models, where holon 
excitations may be stable. But also the cases of very 
weak and very strong interaction allow long-lived holon 
excitations. In this section, we shall therefore develop the 
theory for holon edge singularities assuming the holon is 
stable. 

Let us denote the holon threshold in the spectral func- 
tion by e c {k). By definition, the configuration giving rise 
to this threshold contains a spinon at the Fermi point as 
well as a holon excitation of energy e c (k). A first step 
towards the calculation of the corresponding threshold 
exponents is generalizing Sec. [IV] to calculate the phase 
shifts caused by such a holon impurity. The derivation 
is very similar to the spinless case^i The Hamiltonians 
Hq and Hd describing, respectively, the LLs at the Fermi 
points and the holon impurity arc given by 

K c {ve c ) 2 + -L(V0 C ) 2 

Ac 



H = £ld X 



(V0 S 



dx dl(x)[e c (k) - iv d V]d c {x), 



(69) 



where Vd = de c (k)/dk. Interactions between the holon 
impurity d c and the spinons at the Fermi points do 
not lead to phase shifts because of S'[/(2)-symmetry: a 
density-density interaction of the type p as (x)d\{x)d c (x) , 
being linear in spinon density would violate spin- 
up/spin-down symmetry. Therefore, the interaction term 
only contains holon-holon interaction and becomes in 
bosonized form 



Hi, 



dx 
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d\d c . (70) 



Removing the interaction term by means of a unitary 
transformation as in the spinon case leads to phase shifts. 
In order to distinguish them from the spinon phase shifts, 
we label them x±- In analogy to Eq. ([4*0]) . these phase 
shifts are determined by the equations 

(V[ T V'r)KT 12 = -X-(v d + v c ) ± X +(vd - v e ). (71) 

In order to express the interaction potentials V L R in 
terms of measurable quantities, we consider again the 
variations of the energy of the system with respect to a 
variation of the density and to a Galilean boost using the 
Hamiltonian Hq + Hd + Hi nt . A uniform density varia- 
tion Sp c leads to a nonzero (V0 C ) = —ir8p c /\/2 and thus 
shifts H int by 
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vk 



2V2 



(72) 



The density variation also causes a change in the holon 
Fermi momentum but this does not lead to an energy 
shift in Hq and Hd- Calculating the same shift in energy 
using the definition of e c (k) [see Eq. ([4i~jl ] and equating 
both expressions leads to 



Vk + K _ de c (k) irv 



2V2 



dp c 



w- (73) 



Next, we consider again the energy shift due to a uniform 
change of momentum. For a system moving at velocity 
Su, this leads to a finite (V0 C ) = \/2m5u. Equating the 
corresponding shift in Hi nt with the shift calculated from 
the definition of e c (k) [see Eq. (|45p] leads to the second 
relation, 
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(74) 



The phase shifts can now be determined from Eq. (|71l) . 
The result is 
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(75) 



This result is therefore a direct generalization of the re- 
lation derived for the spinless case in Ref. [2l| . The sym- 
metry of the holon edge e c (k) — e c (— k) leads to the sym- 
metry x±(fc) = — x=p(— k). Note that in contrast to the 
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spinon momentum, the holon momentum is not bounded. 
Therefore, the relations ([75} hold for arbitrary k. 

The phase shifts (|75j) reproduce the correct universal 
phase shifts in the vicinity of k — > fc^. Close to fcp 5 one 
can expand e c (k) = v c (k — kp) + [k — kp) 2 / (2m*), where 
to* is an effective mass which will generally be different 
from the bare mass m of the physical particles. Using 
this form of e c (k) it may easily be checked that 



X-{k — > kp) = 5- 



(76) 



where (S_ c is defined in Eq. fj 16[) . In order to verify the 
correct behavior for x+, it is necessary to determine the 
effective mass to*. This can be done by bosonizing the 
Hamiltonian with quadratic spectrum and considering 
the response of the system to a density variation. The 
result is similar to the spinless case^ 5 - and reads as follows 
for a Galilean-invariant system, 



1 _ VZKcdvc 

TO* 2?7 



(77) 



VII. LIMITING CASES 

A. Strongly interacting fermions 

For strong repulsive interactions, it becomes increas- 
ingly difficult for fermions in one dimension to pass each 
other. The excitations can be separated into charge and 
spin parts and the spin part can be modeled using a 
Heisenbcrg Hamiltonian, H s — J~^2j Sj ■ <Sj+i, where Sj 
denotes the spin density on site j of the lattice. A strong 
finite-range interaction leads to an exponentia l 58 ' 66 " — 
suppression of J and, as a consequence, the spinon spec- 
trum becomes almost flat, |e s (fc)| J$ J. In this case, the 
phenomenological phase shifts defined in Eq. (|48| are 



A<5 ±c (fc) 
2tt 



2V2 k F 



(80) 



The edge exponents ± in the range (2n — l)kp < k < 
(2n + l)kp can then be determined from Eqs. (|55|) and 
(|56[) using the proper selection rules for m±. They are 
given by 



Using this effective mass, it can be shown that x+(k —> 
kp) = S +c . 

Power-law singularities of the spectral function may 
appear at e c (fc) as well as at the shifted holon lines 
e c {k — 2nkp) for n G Z. The configurations giving rise to 
singularities at uj w e c {k — 2nkp) contain a holon which 
carries almost the entire energy ui of the incoming parti- 
cle (or hole), as well as \n\ additional low-energy particle- 
hole pairs with total momentum 2nkp. The calculation 
of these edge exponents is analogous to the calculation in 
Sec. IV Bl The fermion annihilation operator is projected 
as 



ipt = e lkx d c F Rc ip Rs F R3 x <f>2nk f 



(78) 



where $2nt F is given by Eq. ([53]) . The edge exponents 
can now be found by bosonizing <&2nk F and ij) Rs using 
Eqs. © and (TTlD. 

Assuming the holon energy is well-defined, the spectral 
function displays a power-law behavior near the holon 
spectrum, A(k,u>) cx [a; — e c (k — 2nfcj^)] _A V- with 
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(79) 



Note that m_ has to satisfy the selection rule in Eq. (f5T|) , 
so the leading exponent is always reached for to_ (to_ + 
1) = 0. The edge positions and the labels of the ex- 
ponents in the different regions of the (fc, w)-plane are 
illustrated in Fig. 
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For k — > kp, this reproduces the universal exponents (|2ip 
of the LL theory. 

In general, the holon branch of excitations is broadened 
by possible decay via the creation of pairs of spinons. 
However, in the limit J — > 0, we expect the interactions 
between the holon and spinon branches to vanish^ In 
this limit Eqs. (|75jl and (|79|) are applicable. The result- 
ing exponents depend on details of the interaction which 
ultimately defines the dependence of the holon spectrum 
on momentum and density. Below we consider the special 
case of a Hubbard model, in which the existence of holon 
mode is protected by integrability, and the interaction 
becomes strong at low electron filling factors. 

The Hubbard model is described by the Hamiltonian 

^Hubbard = ~t^2 [lpl a 1p( n+ l)a + h.C.] 

7 1(7 

+(u/2)^UMi^m, (82) 

n 

where ip na annihilates a fermion of spin a =^, ^ at the nth 
lattice site, t is the hopping amplitude between neigh- 
boring sites and U is the on-site interaction strength. 
For large interaction, double occupancy of a single lat- 
tice site is suppressed and the Hubbard Hamiltonian can 
be mapped onto a t-J model,— where the J-term de- 
scribes a coupling between neighboring spins with an ex- 
change coupling J = At 2 /U . For large U, the spinon 
spectrum e s (k) collapses. Moreover, U — > 00 ensures 
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that double occupancy of a single lattice site is forbid- 
den and the charge degrees of freedom of the spinful in- 
teracting fcrmions behave largely like spinless nonintcr- 
acting fcrmions of the doubled density. Therefore, the 
charge part can be regarded as noninteracting particles 
with Fermi momentum 2kp and hence Fermi velocity 
V c = 2vp. Going back to the continuum case, this cor- 
responds to a Luttingcr parameter K c = vp/v c = 1/2. 
The edge exponents become 



even,± 



1 l / k 

T 2 " 8 \fo 



(83) 



They coincide with known results determined from the 
finite-size corrections^ of the exact solution of the Hub- 
bard model using the Bethe ansatz and from the fac- 
torization of the exact solution into charge and spin 
parts' - — 

Being an integrable model, it is expected that the Hub- 
bard model has stable holon excitations and that the 
power-laws in the spectral function at the holon mass 
shell are resolved. The exponents are given by Eq. ([79| . 
were the phase shifts x±(k) have to be determined us- 
ing the holon spectrum e c (k) for the strongly interacting 
Hubbard model. From Eq. (|77[) . one finds that the ef- 
fective mass for K c = 1/2 is equal to the bare mass, 
to* = to. Moreover, for U — > oo, the charge sector be- 
haves like noninteracting spinless fcrmions with Fermi 
momentum 2kp. Hence, the holon energy as a func- 
tion of the holon momentum k c is given by ui c (k c ) = 
£;;?/(2TO)-(2fc F ) 2 /(2m), see Fig.[S This is also confirmed 
in the Bethe ansatz solution. For k > 0, the edge position 
in the spectral function is given by e c (k) = uj c (k + kp). 
The symmetry e c (k) = e c (|fc|) then fixes the thresh- 
old for arbitrary k. This threshold position has been 
found numerically in (75| . Using Eq. (|75j) . one readily 
finds X±/(27r) = 1/4 independent of k. Therefore, for 
io ~ £ c (k), we find A(k,u>) cx [uj — e c (fc)] _A1 ° - with 

M§,- = | (84) 

In the limit k — > kp, this reproduces the universal result 
Eq. (|25|) . A schematic picture of the spectral function 
for the strongly interacting Hubbard model is shown in 
Fig.lU 



B. Weakly interacting fermions 

For noninteracting fcrmions with dispersion e(fc), 
A{k,ui) = 5[ui — e(k)]. As we have seen in the previ- 
ous sections, for |w| > |e(fc)| interactions generally turn 
the singularity at the mass shell into a power-law. More- 
over, for > fci?, the edge of support no longer coin- 
cides with the mass shell but, in the limit of vanishingly 
small interaction, with the shifted and inverted threshold 
both = —e(k ± 2kp), see Fig. [TT] In this section, we shall 
calculate the spectral function using perturbation theory 
in the interaction. 
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FIG. 9: (Color online) Position of power-law singularities in 
the spectral function A(k, lj) for the strongly interacting Hub- 
bard model. At the main holon branch e c (k), power-law sin- 
gularities have an exponent pf) t _ regardless of the sign of uj. 
Weaker power-laws with exponents ^i^- can be found at the 
shifted holon mass shell e c (k + 2nkp) for n £ Z. 



For this purpose, let us start from the definition of 
the spectral function in terms of the retarded Green's 
function, A(k,uj) = -\ Im G rot (fc, uj). Due to SU(2)- 
symmetry, all correlation functions are independent of 
the spin orientation, so the spin index was dropped. The 
exact retarded Green's function of the interacting system 
can be expressed in terms of the retarded self-energy as 



G ret (/c,w) 



1 



e(fc)-E ret (/c,w) 



(85) 



The self-energy 57 ct (fc,aj) will be calculated pcrturba- 
tively using the Hamiltonian H = Hq + Hi nt , where 

ffo=EE#l(*)w*). ( 86 ) 

k CT=f,4. 

^Uki)^Uh)Vi(k 3 )Mk2 + h)Mh - h) 
-V4(fci)4(fc2)^ 2 (fc 3 )^(fc 2 + h)^ T (h - k 3 ) 

The free fcrmion spectrum is c(k) = ( k 2 — kp)/ (2m) . and 
ip„ (tpl) denotes the annihilation (creation) operator for 
a physical fcrmion of spin a =j" , J.. These operators obey 
the fermionic commutation relations, {ip a (k),ipl(k')} = 
Sa-rSkk'- For system length L, the momentum k = 2im/L 
(n G Z) is quantized due to the periodic boundary condi- 
tions, but we shall take the limit L — > oo in the following. 
The term Hi nt is the most general two-particle interac- 
tion term allowed by SU (2)-symmctry and translation 
invariance. In terms of the charge density p c (x) and the 
spin density S(x) defined in Sec. IV C\ one can write it as 



(87) 



H int = - J dxdy[U p (x - y)p c (x)p c (y) 
+ U a (x-y)S(x)-S(y)]. 



The interaction potentials of charge and spin densities 
are related to Vi. 2 by Vi = U p - U a /4 and V 2 = -U a /2. 
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The sclf-cncrgy will be calculated pcrturbatively in the 
interaction strengths Vi <2 - Separating real and imaginary 
parts of the self-energy, the spectral function reads 



A(k,L0) 



1 



Im E r 



7T [w - e(fc) - Re £ rct ] 2 + [Im £ rct ] 2 ' 



(88) 



The first-order contribution to £ rot (fc,w) can be ac- 
cessed most conveniently by calculating the first-order 
self-energy H(k,iu) n ) in imaginary time and then per- 
forming an analytic continuation iuj n — > uj + iS in order to 
translate this to the retarded self-energy. The only con- 
tributions in the first order are the well-known Hartree 
and Fock terms. The result is energy-independent and 
reads 



£ ret «(fc) 



dq 



ik-k F 27r 
2V 2 (q) + V 2 (0) 



[V 1 (q)-2V 1 (0) 



(89) 



Because the first-order term is real, it merely leads to a 
shift of the edge position. As the energy uj is measured 
with respect to the chemical potential, the new edge po- 
sition is shifted to 



e'O) = e(k) - £ ret «(fc) + E ret W(fc F ), 



(90) 



and the spectral function up to first order in the in- 
teraction strength remains a (5-function, A^(k,u>) = 
6[w-(?(k)]. 

The second-order self-energy could in principle also be 
calculated using the imaginary-time Green's functions. 
In view of the calculation of the power-laws in the spec- 
tral function, however, our primary interest is in its imag- 
inary part. The most convenient way to calculate the 
imaginary part of the self-energy to second order is pro- 
vided by Fermi's Golden Rule, 

-ImE ret < 2 )(fc, W ) = (f\ R ™t l*> l 2 ^ - e f ). 

(91) 

The initial state \i) has energy and momentum uj and 
k, respectively. The sum over final states |/) is over a 
complete basis of the Hilbert space and e/ denotes the 
energy of the state |/). The (5-function reflects energy 
conservation. 

Let us investigate first the particle sector and as- 
sume k > hp. In this case, the initial state is given 
by \i) = Vi(fc) | ^-S") , where \FS) denotes the Fermi 
sea of spin-up and spin-down particles where all states 
with |fc| < kp are filled. The most general final state 
yielding a nonzero matrix element is given by |/) = 
ip a ,(k — q)4>l,{p + q)ip T (p) \FS), so Eq. (f9"Tj) measures the 
decay probability of a single particle with momentum 
k and energy uj (not necessarily on mass shell) via the 
creation of a particle-hole pair with momentum q. The 
physical process is depicted in Fig. [TUJ 

Let us first calculate (|9"T|) for k > kp for energies 
slightly above the mass shell, ui ^ e(k). In this case, 




FIG. 10: Feynman diagram for the imaginary part of the 
second order self-energy. The straight lines denote fermions 
with momentum and spin quantum numbers. The wiggly line 
depicts the interaction with momentum exchange q. 



energy and momentum conservation imply that either 
q k + kp or q J$ 0. Both processes are identical since 
they can be mapped onto each other by exchanging the 
momenta p + q and k — q of the outgoing particles. Phys- 
ically, for the quadratic spectrum e(k), a right- moving 
particle with k > kp and ui > e(k) can only decay by cre- 
ating a particle- hole pair near the opposite Fermi point. 
For small uj — e(k), 



ImE rc '( 2 >()t, W ) 



2U* s (k + k F ) 



6[uj~e(k)}- 



e(k) 



1 (v + Vf) 2 ' 
(92) 

where we used v = k/m and vf = kp/m. For interaction 
potentials V\^ 2 {q) which do not vary appreciably on the 
scale Aq = m[us — e(k)]/(k + kp) near q = k + kp and 
q = 0, the effective interaction vertex is given by 



[Vi(0)- 



'V 2 (k)} 2 + [V 2 (0) + V 1 (k)} 2 
V 2 (k)][V 2 (0) + V 1 (k)]. 



(93) 



Next, let us focus on k > kp but energies slightly be- 
low the mass shell, uj ^ In this case, energy and 
momentum conservation are satisfied for q J$ k — kp and 
q ^ 0. As previously, both types of processes are related 
by exchanging the momenta on the outgoing lines. The 
result reads 

- In^Hk, uj) = 2 ^' fcf V ^)-,]^|. 

7T [V — VfY 

(94) 

The spectral function near the mass shell can now be 
obtained from Eq. (|88|) . Due to the analytic struc- 
ture of the retarded Green's function, the real part can 
be estimated from Kramers-Kronig relations. Close to 
the mass shell, these predict a logarithmic divergence 
Rc£ rct ( 2 '(fc, uj) cx U 2 s [uj ~ e(k)} ln{[w - e(k)]/e F } where 
e F is a high-energy cutoff of the order of the Fermi energy. 
This indicates an expected breakdown of the perturba- 
tion theory in a narrow vicinity of the spectrum. For 
a sufficiently weak interaction and at fixed (k — kp) 2 /m 



there arc domains of energy uj where Eqs. (|93p and 
are valid, while Re £ ret ( 2 ) can be dispensed with in com- 
parison with \uj — e(k)\. In these domains, the spectral 
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FIG. 11: (Color online) Singular lines in the spectral func- 
tion A(k, uj) for a weakly interacting systems with spectrum 
e(k) — (k 2 — kp)/(2m). Up to second order in the interaction 
strength, singularities appear at the mass shell to « e(k) for 
all k as well as at to ~ —e{k =p 2/cf) for k ^ ±&f- In the 
configurations giving rise to the singularities, the energy is 
carried by a spinon along the red (solid) line or by a holon 
along the green (dashed) line. 



function reads 

A(k,uj) 



2U 2 B {k + k F ) 0[ 
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tt 2 (v + vf) 2 oj — e(k) 
2U 2 s (k - k F ) 0[e(k) - lu] 
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(95) 



Therefore, comparing this to the general structure 
A(k,iij) oc [lu — e(fc)] _Ato -, perturbation theory predicts 
an exponent /i§ _ = 1 in the vicinity of the mass shell. 
The singular lines of the spectral function in the weakly 
interacting case are depicted in Fig. [TT] 

Next, we calculate the spectral function for kp < k < 
3kp near the edge of support, uj w Lu t h = —e{k — 2kp). 
It turns out that for small lu — ujthi due to energy and 
momentum conservation, the only allowed momentum 
exchange must satisfy q ;$ k — kp. To lowest order in 
lu — Uthi one finds 



■ImE ret ( 2 '(fc,i 



[UJ - Uth\ 



uj - ujth 



7T ' (v — Vf) 2 

with the effective interaction vertex, 

V c % = - k F ) + V 2 (k- k F )] 2 . 



(96) 



(97) 



The same caveats as previously about the applicability 
of perturbation theory apply also near the threshold. In 
the range of applicability, one finds the spectral function 
from Eq. 



A(k,u) 
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"off 
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(lu - Lu t h)(u - wth)- (98) 



Therefore, the threshold exponent predicted by pertur- 
bation theory is [i\ + = — 1. 

The second-order self-energy can be calculated for 
— kF < k < kp in a similar fashion. In this case, the 
initial state reads \i) = ip a (k) \FS) while the final state is 
given by |/) = V4(p)Vv(j>+<7)Vv(fc-<7) \FS). The calcu- 
lation of the matrix elements (/| Hint |*) is performed as 
previously but now leads to different integration ranges 
for the momenta p and q. For uj w e(k) < 0, a hole 
with momentum k can relax by creating a low-energy 



particle-hole pair near either the left or right Fermi point 
by transferring the momentum k — kp or k + kp, respec- 
tively. For |fc| < fci?, the spectral function near the mass 
shell reads 



A(k,oj) 



U 2 s (k + k F ) U 2 s {k-k F ) 



tt 2 (v + vf) 2 tt 2 (v — vf) 2 



0[e(k) 
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(99) 



This corresponds to an exponent _ = 1. For |fc| < kp, 
the mass shell in the limit of vanishing interaction co- 
incides with the edge of support. However, note that 
second-order perturbation theory is insufficient to explain 
the finite values of A(\k\ < kp,uj) near the shadow band 
u> w — e(fc) > 0. The final state in this region contains 
an additional particle- hole pair with momentum ±2kp 
and is thus not captured by the second-order calcula- 
tion. Note also that unlike the case of spinless fermions, 
here U c s(0) ^ 0, indicating problems with the perturba- 
tion theory in the vicinity of Fermi points; we have to 
require \k — kp\/m U c g. We note that the difference 
between the holon and spinon velocities, v c — v s oc U c s, 
so the latter requirement may also be viewed as the con- 
dition that the spectrum curvature term in the fermion 
energy is appreciable, (fc — kp) 2 jra ^> (v c — v s )\k — kp\. 
Because the particle-hole pair in the final state has to 
respect this constraint, perturbation theory can only be 
used to evaluate A(k, u>) for energies uj away from the true 
singularities, \lu — e(k)\, \lu — uu t h\ 3> {v c — v s )\k — kp\. 

In order to calculate the exponents in closer vicinity 
of the thresholds, we use the phenomenological relations, 
Eqs. P5]) and (|75|) . First, let us explain how to decom- 
pose a physical fermion into spinons and holons in the 
weakly interacting limit. Wc use the Bcthc ansatz solu- 
tion of the Yang-Gaudin model^i^ in the limit of zero 
interaction and compare the thresholds e(k) and uu t h to 
its spinon and holon spectra.— For k < kp, the singular- 
ity at uu « e(k) is created by configurations which contain 
a holon at its Fermi point while the energy e(k) is carried 
by a spinon on mass shell. For k > kp, the excitations at 
lu w Lu t h arc spinon excitations in the same sense. Hence, 
the general observation that for repulsive interactions the 
edge of support of the spectral function corresponds to 
spinon excitations continues to hold in the noninteracting 
limit. In contrast, for k > kp, the threshold atw« e(k) 
corresponds to a holon excitation with energy e(k), while 
the spinon rests at its Fermi point. The singular lines of 
A(k, lu) in different regions of the (fc,a;)-plane along with 
the nature of the respective excitations are depicted in 
Fig.HU 

Away from kp, the physical spectrum e(k) is only 
weakly affected by the interactions, so we can calcu- 
late the exponents in this region using the noninteract- 
ing spinon and holon spectra. For < k < 2/cf, ex- 
cluding again a domain \k — kp\ ~ mU c ff(0) around 
the Fermi points, we can use e s (k ^ kp) = vp(k — 
kp) T (k — kp) 2 /(2m) and calculate the exponents from 
the phenomenological relations Eq. For weak in- 

teractions, the Luttingcr parameter can be expanded 
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as K c = 1 — 5K C , where 6K C —> + . Everywhere at 
|ft| < kp, except the narrow vicinities of the Fermi points, 
| ft ± kp\ ~ mU e g, the threshold exponent is given by 
lUg _ ps 1 — 0(8 K 2 ) and thus is compatible with the per- 
turbative result (|99|) . For ft > kp and uj ps uth, the 
phenomenology yields the exponent /xf + = 0(6 K c ) at 
the threshold, different from the exponent Eq. valid 
away from the threshold. 

Similarly, the holon spectrum at \k — kp\ 3> mU e fi can 
also be approximated by its nonintcracting limit, e c (ft > 
kp) = vp(k — kp) + (k — kp) 2 1 (2m). A calculation of the 
phenomenological phase shifts and the exponent using 
Eqs. (Z5J) and (75) then leads to j/g _ = 1/2 - 0(6K 2 ). 
Note that this coincides with the exponent at the holon 
mass shell predicted for k m kp by the linear LL theory, 
see Eq. ([21]) . But it is different from the perturbativc 
exponent /xg _ = 1 in Eq. (|95[) which is valid away from 
the mass shell. 

We conclude that the lowest-order perturbation the- 
ory performed here is only able to predict the behavior 
of A(k, u) away from the true edge. Closer to the edge, 
the nonperturbative spin-charge separation becomes im- 
portant and the correct exponents can be derived using 
the phenomenological relations. 



| ft — ftp |, possibly contradicting to the notion of a well- 
defined holon branch at k — > kp. However, as we show 
below, the decay rate of a holon with small (measured 
from the Fermi point) momentum must scale to zero with 
k — > kp faster than \k — kp\ 3 . 

In order to elucidate the possible decay processes for 
holons, it is convenient to start again from a descrip- 
tion in terms of refermionized quasiparticles. The band 
curvature of the physical fermions leads to interactions 
between the quasiparticles. Away from the Fermi points, 
it is advantageous to classify the interaction processes by 
their relevance in the RG sense and to consider all pos- 
sible interaction operators which arc allowed by SU(2)- 
symmetry and Galilean invariance. Due to its built-in 
S'[/(2)-symmetry, non-Abelian bosonization^ is a conve- 
nient tool to achieve this. Expressed using the left- and 
right- moving holon densities J a (x) and spinon densities 
J a {x) (a = L,R), the Hamiltonian © of the linear LL 
reads H = H c + H s , where 

H c = 2nv c J dx[J 2 R (x) + 



H„ = 



27TU, 



dx[.P R (x) + Jl(x)} 



n, 



(100) 



VIII. HOLON RELAXATION 

The spinon and holon excitations of the linear LL the- 
ory are eigenstates of the Hamiltonian Hq [see Eq. ([6])] 
and are thus predicted to be stable. As a consequence 
of the infinite spinon and holon lifetimes, the singulari- 
ties of the spectral function and other dynamic response 
functions at the spinon and holon mass shells are char- 
acterized by true power laws. 

In the case of nonlinear spectrum, bosonization leads 
to spinon-spinon and holon-holon interactions, as well 
as to spinon-holon interactions. This may give rise to 
finite lifetimes for holons (we assume that spinons remain 
the lowest-energy excitations). In the spectral function, 
such a finite lifetime generally leads to smearing of the 
singularities^ and is thus directly measurable. 

For spinless fermions, the effects of a nonlinear spec- 
trum on the particle and hole lifetimes were analyzed 
in [20l | . For weakly interacting fermions with a generic 
short-ranged repulsive interaction potential V(k), it was 
shown that the decay rate of a particle with momentum 
k ~ kp scales as T oc (k — kp) s . The exponent here comes 
from a limitation on the phase space available for the de- 
cay combined with the q 2 momentum dependence of the 
effective interaction at small momentum transfer q. An 
important exception is the case of integrablc models,— ~— 
where these calculations predict infinite lifetime for all 
momenta. Away from the limit of weak interactions, the 
rate is limited by the phase space argument only, yielding 
the decay rate scaling^ T oc (k — kp) 4 . 

A similar phase space argument applied to the decay 
of a holon with creation of two spinons would lead to T oc 



The operators J a (x) are related to the physical charge 
density by p c (x) = 2\/K~ c (Jl(x) + Jr(x)). This Hamil- 
tonian emerges at the low-energy RG fixed point and is 
valid in the narrow-band limit. The leading correction 
for increased bandwidth is an interaction between left- 
moving and right-moving spin densities^ 



H g = —2TTv s g / dx Jr(x) ■ Jl(x). 



(101) 



Note that when expressed in terms of the Abclian spinon 
fields (f> s and 8 S , the operator H g generates the sine- 
Gordon term ([7]). The band curvature of the physical 
fermions leads to interaction operators which are cubic 
in spin and charge densities, 



H„ 



H K = 



4tt2 

ir 

4tt2 
3 



J dx [J7_(4 + J£) - r 1+ (J 2 J L + J 2 J, 



dx 



k-(JrJ r + JlJ[) 



He 



4tt 2 C 



dx (J L + Jr)Jr ■ J L - 



(102) 



Note that these operators represent all cubic terms which 
are compatible with Sf7(2)-symmetry. In particular, this 
symmetry prohibits terms linear in the vector opera- 
tors J a (x). Interaction operators containing quartic and 
higher-order terms in J a (x) and J a (x) do exist but their 
contribution is sublcading for small bandwidths. 

The prefactors g, £, k± and r)± can be fixed phe- 
nomenologically by relating them to other observable 
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quantities. For this purpose, let us investigate the varia- 
tion of the interaction operators as a response to a uni- 
form variation of the charge density p c {x) — > p c (x) +6p c . 
This variation shifts J a {x) — > J a (x) + 5J, where 5 J = 
Spc/^y/Wc), and thus creates the following additions to 
the Hamiltonian, 



SH K 



SJ 



(3 V -- V+ ) / dx{J R + J 2 L ) 



4r/ + / dx J L J R j, 
: -^-6J{k_ + k+) I dx{Jl + ,P L ), 



{28 J) / dx J r -Jl- 



(103) 



Combined with H c , the Hamiltonian SH V leads to a 
renormalization of v c and K C J& The shifts SH K and 
5H^ modify H s and H gi respectively, and thus renor- 
malize the parameters v s and v s g. Expressing the den- 
sity variation as a variation of chemical potential using 
dpc/dp = 2K C / (irv c )r one hnds^ 



(104) 
(105) 



In a similar way, the difference k_ — n + can be related 
to the mass m of the physical fermions by considering a 
charge current variation of the Galilean-invariant system. 
One finds^^o 




K- 



(106) 



It is known that upon a bandwidth reduction g flows 
logarithmically to zero^ Assuming the initial bandwidth 
to be of order kp, for a smaller bandwidth of order 
k the effective coupling constant will flow to g(k) = 
1/ \n[kp/(k — kp)]. As the chemical potential p is pro- 
portional to kp, the derivative dg/dp oc ~g 2 jkp. The 
derivative dv s /dp, on the other hand, remains finite for 
small bandwidths. Therefore, in leading logarithmic ap- 
proximation, dg/dp can be neglected and the coupling 
constants k± and £ can be related as 



c 



dv s 3 
^p- = -2 9{K - 



(107) 



For repulsive interactions, the excitation of lowest en- 
ergy for a given momentum k is a spinon with energy 
e s (k). Due to the shape of the spinon spectrum (see 
Fig. [7]), the absolute value of the velocity de s (k)/dk 
reaches its maximum v s near the Fermi points. As a 
consequence, spinon relaxation by creation of low-energy 
spinons is ruled out by energy and momentum conser- 
vation. Similarly, because v c > v s > de s {k) / dk, decay 



of spinons by the creation of holons is also forbidden. 
Therefore, spinon excitations are stable. 

Holons, on the other hand, can relax via the creation 
of low-energy spinons. Let us investigate the decay of 
an initial state |i) = \k) c |0) s which contains an addi- 
tional holon with momentum above the Fermi edge and 
no spinon excitations. Relaxation of the holon to a mo- 
mentum k' < k can happen via the creation of two spinon 
density excitations with momenta < and qpi > 0. 
This final state will be labeled |/) = \k') c \qL,QR) s - For 
momenta k close to the Fermi point, momentum and en- 
ergy conversation for this process read 



k = k' + q R + q L , 
v c k = v c k' + v s (q R -q L ), 



(108) 



and have nontrivial solutions (k ^ k') for v c > v s . 

The holon lifetimes associated with this decay channel 
can be calculated using Fermi's Golden Rule. Two com- 
binations of operators from the interaction terms (|102[) 
have a nonzero matrix element between the states |i) and 
|/). To first order in the interaction, (f\H^\i) is the 
only such term. To second order, only (/| H g H K \i) and 
(/| H K H g \i) are nonzero. 

The calculation of these matrix elements is greatly 
simplified by the fact that the initial state \i) contains 
no excitation in the spin sector and thus corresponds to 
an 5[/(2)-invariant singlet state. Therefore, one can use 
J a (x) ■ J p {x) |0) s = 3J^(x)J^x) \0) s for a,P € {L,R}. 
Moreover, all interaction Hamiltonians H g , H rp H K and 

H{ conserve the total spin S = J dx[J R {x) + Jl(x)]. 
This can be verified by calculating the commutators 
using the 5f/(2)-Kac- Moody algebra for the operators 
J a (x)£^ Therefore, the spin sector remains in a singlet 
state even if acted on by the interaction Hamiltonians. It 
means that only the z-componcnts of the spin operators, 
Ja(x), are needed for the calculation of the holon life- 
time. A normal-mode expansion allows one to represent 
the Fourier transform (p) in terms of bosonic operators 
b„ 



J Z M 



47T 



[6(ap)bl + 6(-ap)b. 



(109) 



where L is the system length and a = R, L = +, — . The 
operators b q satisfy the commutation relations [b q ,b* ql ] = 
8qq' . Using the creation operators b q , the spin part of the 
final state |/) can be written as \qL,qn) s = b\ L b qR |0) s . 
For the calculation of the holon lifetime, the following 
spinon matrix elements are needed, 



(q\J^p)\0) s = X l^0(aq)6 p , q . 



(110) 



For the expectation values of holon operators, on the 
other hand, it is convenient to retain the description 
in terms of fermionic quasiparticles where J a (x) = 
\p ac {x) = i^ ac (x)iljac(x) ■ The charge part of the initial 
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and final states can be expressed as \k) c = 4> Rc {k) \FS) C 

and \k') c = ip Rc {k') \FS) C , respectively, where \FS) C de- 
notes the Fermi sea of holons. Then, one finds the matrix 
element, 



z{k'\ Jr(p) \k) c = -Sk',p+k- 



(111) 



The first-order matrix element TV = (f\ \i) can be 
calculated by using these matrix elements. One finds 



L 2 ^ 



(Ql\J Z L ( P )\0) s s (q R \J Z R (p)\0) f 



Xc(k'\J R (-p-p')\k) c 



2^h-k' -q L -q R V\lLqR\ 



(112) 



The matrix elements (/| H g \i) and (/| H v \i) vanish be- 
cause H g and H^ do not couple spinons and holons. The 
remaining first-order matrix element (f\H K \i) — be- 
cause it contains only terms of the form J 2 (x), which do 
not create spinons on opposite branches. 

To the second order, cross-terms of the operators H g 
and H K may couple the same initial and final states as 
above. Since these are second-order terms, the matrix el- 
ements can be calculated using the S'-matrix expansion^ 



T Kg = (f\H K 



E-Hn 



-H„ U) 



T BK = (J\H l 



-H K \i) 



9 E - Hq ' 



(113) 



where Hq = H c + H s is the noninteracting Hamiltonian. 
The energy E denotes the energy of the initial state, 
E = v c (k — hp). As Fermi's Golden Rule imposes energy 
conservation, it ultimately becomes equal to the energy 
of the final state. After Fourier transforming H g and 
H K , going over to bosonic operators using Eq. (|109j) and 
using the spinon and holon matrix elements (|110[) - (|llll) . 
one finds that 



rp _ Sng 
T gK = 0. 



K+)6, 



k-k' - qL -q R \f\qLqR 



(114) 



Other second-order terms exist but they contain higher 
powers of q^ and q R and are therefore sublcading com- 
pared to T Kg for holon momenta k near the Fermi points. 
According to Fermi's Golden Rule the lifetime is 



I/) 



(115) 



where e$ and ej are the energies of the initial state \i) 
and the final state |/), respectively. The sum over all 
final states |/) translates to a summation over the mo- 
menta < 0, q R > and k! £ [kp,k]. It can be 
seen from Eqs. (|112[) and (|114l) that each of the decay 
channels taken individually would lead to a decay rate 



r oc (k — kp) 3 . However, Fermi's Golden Rule (|115j) con- 
tains the square of the sum of the probability amplitudes 
and T Kg . The prefactors of both amplitudes are re- 
lated according to Eq. (|10T[) and one finds Tf + T Rg = 0. 
Therefore, the decay rate vanishes up to terms propor- 
tional to g 2 (k — kp) 3 , in the calculation of T performed in 
the second order— of g = 1/ \n[kp / (k— kp)] . Retaining in 
Eq. (|105[) the derivative dg/d/j, oc g 2 /ep exceeds the ac- 
curacy of our calculation. It is not clear if the evaluation 
of r to order g A would yield zero. Possibly, in that or- 
der the distinction between integrable and non-integrable 
systems emerges. 

In the limit of weak backscattering, V(2kp) <C 
V(0) <C vp, the universal logarithmic dependence for 
g(k — kp) is reached only at very low energies, while 
its bare value g oc V(2kp)/vp is applicable as long 
as [V(2kp)/vp] \n[kp/(k — kp)] ^ 1. In that case, 
Eqs. ([TUB]) and (jTT5|) yield 



r oc tp 



v F 



V(2k t 

Vp 



k — kp 



(116) 



This estimate should be viewed as the result of per- 
turbation theory in V(2kp) in the basis of well-defined 
holon and spinon modes with linear spectrum, which 
sets a limit on holon momenta, k — kp < mV(0) (we 
also used v c — v s ~ 1^(0) in the derivation). Curiously, 
the latter estimate for T at the limit of its applicabil- 
ity, k — kp ~ mV^O), matches the relaxation rate of a 
spinful fcrmion evaluated, in the basis of free fermions, 
by the perturbation theory with respect to the entire 



interaction 
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IX. CONCLUSION 

In conclusion, we have investigated spinful one- 
dimensional interacting Fermi systems at zero temper- 
ature beyond the low-energy regime and calculated their 
dynamic response functions for arbitrary momenta near 
the edges of support in the (fc, cj)-plane. Away from the 
Fermi points, the nonlinearity of the fermion spectrum 
becomes important and the physical properties can no 
longer be explained by the linear Luttinger liquid theory. 
In particular, we shed light on the meaning of spin-charge 
separation away from the low-energy limit. 

The Luttinger liquid theory is based on the assump- 
tion of a linear fermionic spectrum. The cigenmodes 
are spin and charge density waves and the theory can 
be formulated in terms of noninteracting bosonic fields. 
Rcfermionizing these fields in Eq. (|11[) . we introduced 
fermionic quasiparticles, spinons and holons, which con- 
stitute a convenient basis even away from the low-energy 
regime. In contrast to the linear LL theory, a curvature 
of the spectrum of the physical fermions leads to inter- 
actions between spinons and holons. For repulsive inter- 
actions and small |fe| — k F , spinons are the lowest-energy 
excitations. By continuity, we expect that for generic 
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repulsive potentials, the spinon spectrum e s (k) remains 
the edge of support for the spectral function for arbitrary 
momenta. 

We found that the spin-charge separation exists also 
for nonlinear spectrum but in a weaker sense than in a 
linear LL. If a particle with arbitrary momentum k and 
energy u ~ e s (k) tunnels into the system, it creates a 
single spinon with energy close to e s (fc), momentum close 
to k, and velocity Vd — de s (k) / dk. In addition, it creates 
low-energy holon excitations with velocity v c > Vd and 
momenta near the Fermi points, but no additional spinon 
excitations. The created spinon separates in space from 
the charge excitations due to its different velocity. This is 
reminiscent of the conventional spin-charge separation in 
linear LLs. However, in contrast to the linear LL theory, 
such a decoupling only survives for energies close to e s (k) . 

The separation in the momenta of the spinon and the 
low-energy excitations allows us to introduce an effec- 
tive mobile impurity model as a tool for the evaluation 
of measurable dynamic response functions, such as the 
spectral function A(k,ui) and the charge and spin den- 
sity structure factors, S(k, oj) and S zz (k 7 uj), respectively. 
In analogy to the Fermi edge problem, the created spinon 
acts as a mobile impurity and causes a shake-up of the 
Fermi seas of holons. The threshold exponents of A(k, oj) 
can therefore be expressed for arbitrary momenta and 
interaction strengths in terms of scattering phase shifts 
A6± c (k), see Eqs. (|55|) and (|56|) . In the vicinities of the 
Fermi points, A5± c depend only on the Luttinger param- 
eter K c , see Eq. (TlT)|) . which leads to universal values of 
the exponents. For arbitrary momenta and Galilean in- 
variant systems, we related the phase shifts to another 



set of measurable properties given by the dependence of 
e s (k) on k and the charge density p c , see Eq. (|4*5]). Us- 
ing similar considerations, we calculated the threshold 
behavior of the spin and charge density structure factors 
S(k, oj) and S zz (k,oj), see Eqs. §7$ and flBgJ). These re- 
sults are summarized in Table HI The general evolution of 
the spectral function with increasing interaction was con- 
sidered by the analysis of the limits of weak and strong 
interaction. In the former one, we found A(k,oj) pcrtur- 
batively, see Sec. I VII B[ In the latter limit, we utilized 
the phenomenological relations to calculate the threshold 
properties of A(k, oj) from the exactly known spectrum, 
see Sec. IVITaI 

Unlike spinons, the holon excitations do decay as a con- 
sequence of the nonlinear spectrum of the fermions. The 
corresponding lifetime is long for holons with momenta 
near the Fermi points, see Sec. I Villi For integrable mod- 
els, the holon spectrum e c (k) is well-defined for arbitrary 
momenta. We determined the exponents of the spec- 
tral function for oj rs e c (k) by deriving phenomenological 
expressions relating the corresponding scattering phase 
shifts to properties of e c (k), see Sec. IVII 
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